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Abstract: We complete the analysis of transport phenomena in p-wave superfluids within 
gauge/gravity duality, using the SU{2) Einstein- Yang-Mills model with backreaction. In 
particular, we analyze the fluctuation modes of helicity zero in addition to the helicity 
one and two modes studied earlier. We compute a further transport coefficient, associated 
to the first normal stress difference, not previously considered in the holographic context. 
In the unbroken phase this is related to a minimally coupled scalar on the gravity side. 
Moreover we find transport phenomena related to the thermoelectric and piezoelectric 
effects, in particular in the direction of the condensate, as well as the flexoelectric effect. 
These are similar to phenomena observed in condensed matter systems. 
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1. Introduction 

Gauge/gravity duality has become a valuable tool for gaining insight into the physics of 
many different strongly coupled theories and, in particular, is being used to successfully 
describe their hydrodynamical behavior, with the prospect of making contact with sys- 
tems found in nature. Recently, many new features within hydrodynamics have been 
discovered using gauge/gravity duality: For instance, the importance of anomalies for rel- 
ativistic hydrodynamics as applied to quark-gluon plasma first appeared in the context of 
gauge/gravity duality [1, 2]. Subsequently, in [3-7] and [8, 9] it has been realized by ther- 
mal field theory computations and from general hydrodynamics arguments that anomalies 
induce modifications in the constitutive relations of relativistic hydrodynamics. Moreover, 
in [10-12] chiral anomalies have been shown to give rise to non-dissipative transport coef- 
ficients. Anisotropy has been included by considering the backreacted holographic p-wave 
superfluid (see below) and by means of a position-dependent theta-term, leading to sev- 
eral interesting effects [13-15]. A common feature of these systems is that the breaking 
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of symmetries brings about a richer structure to the theory, so that new phenomena are 
unveiled. 

A very suitable system to study anisotropic hydrodynamics is the holographic backre- 
acted p-wave superfluid, in which the rotational symmetry is broken by a vector condensate 
which may be interpreted as a vector meson. This system has been studied in [16-21]. It 
involves a finite SU (2) charge density or isospin density. In the present article we present 
the study of the remaining hydrodynamic modes that were not accounted for in [19] and 
describe the corresponding new transport properties. 

In [18, 19], the helicity two and one fluctuations have been analyzed. It has been found 
that the helicity one modes lead to contributions to the viscosity tensor whose ratio with the 
entropy density is non-universal at leading order in the 't Hooft coupling and N . These 
contributions are temperature dependent and satisfy the viscosity bound, rj/s > 1/4-k. 
This is in contrast to the O-teim model of [14] where the usual viscosity bound [22-24] 
is violated [25]. This happens already for Einstein gravity, violations of the bound by 
Gauss-Bonnet terms have been studied in [26]. 

The Einstein- Yang-Mills model may also be used as a starting point to derive similar 
universal relations, such as the holographic realization of Homes' law [27] of condensed 
matter physics. Furthermore when considering flnite SU{2) magnetic fields, the system 
admits more than one possible solution (or state), but similarly to the holographic super- 
fluid at flnite SU{2) density, only one is physically realized, determined by the lowest free 
energy. A magnetic fleld generates an Abrikosov lattice [28] in a superconductor, which 
becomes the preferred state if the magnetic field is sufficiently large. 

The Einstein- Yang-Mills system we consider in this publication is motivated by the 
D3/D7 setup [29], which allows for temperature and matter in the fundamental represen- 
tation to be added to the system. Holographically, the Hawking temperature of a black hole 
geometry coincides with the temperature T of the dual thermal fleld theory. A chemical 
potential // can also be introduced by placing a non-vanishing boundary condition upon 
the bulk gauge fleld. Given these ingredients, it is possible to do thermodynamics, since 
each solution labeled with T/ /i describes a different thermal state of the dual fleld theory. 

In this paper we consider a superfluid generated by a flnite SU{2) density (for an ex- 
tensive study of this background see [16, 19]). In this case, the temperature determines 
the preferred state, i.e. at some critical value Tc the system undergoes a phase transition 
between the normal (T > Tc) and superfluid (T < Tc) states. We will be interested in the 
superfluid phase and study its transport properties. To do so, we consider fluctuations in 
a backreacted holographic p-wave superfluid theory deflned in an AdSs geometry with an 
SU{2) Yang-Mills gauge fleld. The boundary condition that flxes the chemical potential, 

lim At = fi, (1.1) 

breaks explicitly this SU{2) symmetry, leaving a U^l)^ gauge symmetry. Depending on 
the values of the system may present a superfluid or a normal phase, with the order of 
the phase transition being controlled by a parameter a, which on the gravity side measures 
the effect of the gauge flelds onto the geometry. The superfluid state is thermodynamically 
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preferred at low temperatures compared to the chemical potential, and the transition to 
this phase is characterised by the formation of a vector condensate {Jf), as opposed to 
the case of an s-wave superfluid, in which a scalar field condenses. The vector condensate 
designates a particular direction both in momentum and flavor space, and as a consequence 
the spatial rotational SO{'i) symmetry and the U{\)^ symmetry are spontaneously broken. 
Schematically, this process can be represented as 

SU{2) >C/(1)3 — >'L2 , 

Expl.B SSB ,^ 2\ 

50(3) — > SO(2) . 

SSB 

This is an example of spontaneous breaking of continuous symmetries in gauge/gravity 
duality (first achieved in [30]) to construct holographic superfluids or superconductors. 
This technique was initially developed by breaking Abelian symmetries [31, 32] and later 
adapted to p-wave superconductors/superfluids [33] as in the case at hand, giving rise to 
the first string theory embeddings of these constructions [34-36] . 

We present an analysis of the perturbations of the spacetime metric, /i^jy, and of the 
Yang- Mills field, a°, about the Einstein- Yang-Mills model in AdS^. Due to the breaking 
of the spatial rotational symmetry, these fluctuations can be grouped according to their 
transformation behavior under the remaining SO (2) rotational symmetry around the x- 
axis. In this paper, we present the fluctuations which transform as scalars under this 
group. Vector and tensor fluctuations have been studied in [18, 19]. To make the equations 
tractable, we set the spatial momentum A; = 0. This simplification leads to an additional Z2 
symmetry under which the scalar fluctuations can be characterized further. We end up with 
two distinct blocks, the first of which contains, among others, the gauge field fluctuation 

and the metric fluctuation htx, and the second one the diagonal metric fluctuations, a^. 
and a^. 

From the field theory point of view, the corresponding correlation functions are related 
to the thermoelectric effect which correlates charge and heat transport in the direction of 
the condensate, since can be identified with an electric field in the x direction and htx 
with a temperature gradient in the x direction. This effect was studied for holographic 
s-wave superfluids [32, 37] and for p-wave superfluids in the transverse directions [19], but 
to our knowledge, this is the first time this effect has been calculated with backreaction 
and in the direction of the condensate. 

The second block contains, among others, the diagonal metric fiuctuations h^x — hyy and 
the gauge field fluctuations and a^. A field theoretic description of the corresponding 
Green's functions is not fully addressed in this paper and is left for future work. However, 
we know that some of the modes in this block are related to the transport coefficients in 
the viscosity tensor rj^^^K In general, the viscosities of a system are encoded in a rank four 
tensor 77*-''^' which in the most general case has 21 independent components. Due to the 
symmetries of the system at hand, we are left with five independent components of the 
tensor rj^^^'' [38, 39], two of which are shear viscosities, rjxy and r]yz, that were addressed in 
[19]. Two of the remaining components are bulk viscosities and can be set to zero using 
the tracelessness condition for the conformal energy-momentum tensor, leaving one free 
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transport coefficient, denoted by A. While rj and C measure the response of the system to 
deformations due to shear or normal stress, A is related to the normal stress difference that 
is induced by an anisotropic strain. Our holographic computation shows that in the zero 
frequency limit, the ratio of A to the entropy density is finite. Moreover, in the normal 
phase it acquires a fixed value given by 



in the conventions used here. In the broken phase we see a temperature dependence and 
the resulting curve does not fall below the value found in the unbroken phase for any 
backreaction parameter a and for any temperature. 

Since we have completed our analysis of all fluctuation modes in the p-wave system, let 
us now summarize them, as well as the transport phenomena they correspond to: 

1. hyz (helicity two) is related to the shear viscosity r]yz which for all values of T takes 
the universal value r]/s = 1/4tt (see [19]), 

2. hx± (helicity one) is related to the shear viscosity rjx± which shows a temperature 
dependence in the broken phase (see [19]), 

3. The coupling between = it (helicity one) and hx± leads to an effect which 
is similar to the flexoelectric effect known from crystals (see [19]), 

4. Oj^ is related to the "electrical" conductivity a-^-^ (helicity one), and its coupling to 
hf± (helicity one) is related to the so called thermoelectric effect transverse to the 
condensate (see [19]), 

5. <I>4 ~ (helicity zero) is related to the "electrical" conductivity a^^, and its coupling 
to htx (helicity zero) gives the thermoelectric effect in the direction of the condensate 
(see section 4), 

6. <I>3 ~ hxx — hyy (helicity zero) is related to the transport coefficient A found in the 
viscosity tensor ■q^^'''- and its coupling to ^± ~ (helicity zero) shows a behaviour 
similar to the piezoelectric effect (see section 4). 

We see that the study of fluctuations in a backreacted holographic p-wave superfluid 
provides a rich structure of different effects which by using the fluctuation-dissipation 
theorem may be related to well-known transport phenomena in other areas of physics. 

The paper is organized as follows: In section 2, we recapitulate the backreacted holo- 
graphic p-wave superfluid. In section 3, the scalar fluctuations (helicity zero) and the 
corresponding Green's functions on the gravity side are presented. The following section 
4 contains our results for the transport properties of the superfluid and an approach to 
interpreting them from a hydrodynamical point of view. In 5 we present our conclusions. 
Many of the technical details are collected in the Appendices: In A, we discuss the neces- 
sary holographic renormalization, in B we specify the physical fields of the system, in C we 
review the numerical procedure to deal with the coupled equations of motion, and finally 
some general remarks on anisotropic fluids are given in D. 
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2. Holographic Setup and Equilibrium 



The setup used in this paper was aheady described in [16, 19]. Therefore, here we give 
a brief review of its most important properties. We consider SU{2) Einstein- Yang-Mills 
theory in (4 + l)-dimensional asymptotically AdS space. The action is 



S= d^x^ 



^ { n A ^ ^ Z?« TpaMN 



+ S'bdy , (2.1) 



where is the five-dimensional gravitational constant, A = — -p- is the cosmological con- 
stant (with L being the AdS radius), and g is the Yang- Mills coupling constant. It is 
convenient to define 

a = ^ , (2.2) 
9 

which measures the strength of the backreaction. The SU (2) field strength F^^^ is defined 

by 

Fl,^ = dMA% - SnAIi + e'''^A\jA% , (2.3) 



ordinate, and e"^'^ is the totally antisymmetric tensor with e^^^ = +1. The A'^ are the 
components of the matrix- valued Yang-Mills gauge field A = ^^^r^dx*^, where the t"" are 



where capital Latin letter indices run over z,r}, with r being the AdS radial co- 

la 

the SU{2) generators, related to the Pauli matrices by t"" = cr°/2i. Finally, the 5bdy term 
includes boundary terms, namely the Gibbons-Hawking boundary term as well as coun- 
terterms required for the on-shell action to be finite, that will be discussed below. It does 
not affect the equations of motion. 

The Einstein and Yang-Mills equations derived from the above action are 

Rmn + J29MN = t<i {tmn — ■7^Tp^g]\iN^ , (2.4) 

where the Yang-Mills stress-energy tensor Tmn is 

Tmn = ^ (^F^mF'^'^n - ^ffM^vi^^gF'^^^) • (2.6) 

To solve this equations, we use the following ansatze for the gauge field and the metric, 
which can be motivated from symmetry considerations [16, 33] 

A = (j){r)T^dt + w{r)T^dx, (2.7) 

ds^ = -N(r)a(rfdt'^ + -^^dr^ + r^firy^dx^ + r^f(r? (dy^ + dz^) , (2.8) 

N[r) ^ ^ 

where N{r) = — '^"^^'^ + The AdS boundary is at r — t- oo and for our black hole 
solutions we denote the position of the horizon as r^- 

This ansatz is compatible with the well-known AdS Reissner-Nordstrom solution, where 
w{r) = for all values of r. This solution features 'w{r) = 0, so it preserves the 50(3) 
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symmetry and corresponds to the normal phase of the system. There is a second solution 
with non- vanishing w{r), which can only be computed numerically. The second solution 
breaks the rotational 50(3) symmetry and describes the condensed superfluid phase. Due 
to our choice of boundary conditions, this breaking occurs spontaneously. For completeness, 
we state here the coefficients of the expansion at the horizon (in terms of (r/r/^ — 1)"), 

{<^^^o',^o,/o} , (2.9) 

being 0g = in order for A to be well defined as a one- form [40] , and of the expansion at 
the boundary (in terms of (r/j/r)^*^), 

{fi,cPlmlw\jl] . (2.10) 

Note that Wq = 0, otherwise the 50(3) would be broken explicitly instead of spontaneously. 
Besides, we can fix the metric to have asymptotic AdS boundary conditions, so that = 
/q = 1. The fields can be made dimensionless through m(r) — )• r^m(r), (/){r) — )• r/j0(r) and 
w{r) —7- rhw{r), while /(r) and a{r) are already dimensionless. 

In terms of these coefficients we can express the different field theory quantities, such 
as temperature and entropy density, given by 

12-a^^)r;,, s='^rl. (2.11) 




127r 

The field theory expectation values of the dual operators of the different fields are directly 
related to the expansion coefficients. For the charge density and the condensate we have 

{Ji) = -Kd<Pl {JD = K4wl (2.12) 

and for the energy- momentum tensor [41, 42] they are 

^ h h in- \ — h ( b q fb\ In- \ — In- \ — h ( h , , fb 



{ru) = ^ml, {%,) = ^im'^-8f^\ , {Tyy) = {%^) = ^(m', + 4f^ \ . (2.13) 

In [16] it was found that the value of the Yang-Mills coupling constant a determines 
if the phase transition is second order (a < Oc = 0.365) or first order (a > ac = 0.365). 
The quantitative dependence of the critical temperature on the parameter a is given in 
figure 1. The broken phase is thermodynamically preferred in the blue and red regions, 
while in the white region the ground state is the Reissner-Nordstrom black hole. The phase 
transition from the white to the blue region is second order, while the one from the white 
to the red region is first order. The black dot determines the critical point where the order 
of the phase transition changes. In the green region the numerics are unstable. At zero 
temperature, the data may be obtained analytically as described in [21, 43]. 
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Figure 1. This diagram shows the dependence of the order of the phase transition on and Tc/^. 
For a description of the plot see the text. This figure is taken from [19]. 

3. Perturbations about Equilibrium 

In this section we study the response of the holographic p-wave superfluid under small 
perturbations. This analysis is necessary to ultimately compute the transport coefficients 
of the system. On the gravity side, the perturbations are given by fluctuations of the 
metric /iA/Ar(x'^, r) and the gauge field a5^^(x^,r). Thus we are studying a total of 14 
physical modes: 5 coming from the massless graviton in 5 dimensions and 3x3 from the 
massless vectors in five dimensions. Due to time and spatial translation invariance in the 
Minkowski directions, the fiuctuations can be decomposed in a Fourier decomposition as 

/d^fe 

To simplify notations, we drop the hat on the transformed fields which we use from now 
on if not stated otherwise. 

3.1. Characterization of Fluctuations and Gauge Fixing 

In general, we would have to introduce two spatial momenta: one longitudinal to the 
direction of the condensate, /cy, and another one perpendicular to it, k_\_. Thus, = 
(w, fcji , A;^, 0). But introducing a momentum perpendicular to the condensate breaks the 
remaining rotational symmetry 50(2) down the discrete Z2 parity transformation Pj_: 
k±_ — )■ —k_\_, x±_ —7- — x_L. This leads to a mixing of most of the fields making the problem of 
solving the corresponding differential equations unmanageable. Thus we do not study this 
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case further in this paper. However, a momentum exclusively in the direction longitudinal 
to the condensate, or zero spatial momentum, preserves the SO{2) rotational symmetry 
such that we can classify the fluctuations according to their transformation under the 
50(2) symmetry (see table 1). The modes of different helicity decouple from each other. 
The momentum longitudinal to the condensate, however, breaks the longitudinal parity 
invariance Py. In this paper we will set this spatial momentum to zero as well. Therefore, 
we can classify the modes further by there behaviour under the longitudinal parity trans- 
formation P||. Under this transformation the helicity fields are divided into two blocks, 
the first block contains htx, ^1 ^^'^ the second one hu, hxx, hyy + hzz, cix, cl^ 

and Oj. 

In order to obtain the physical modes of the system we have to fix the gauge freedom. 
We choose a gauge where a'^ = and /ijv/r = such that the equations of motion for 
these fields become constraints. These constraints fix the unphysical fluctuations in each 
helicity sector and allow only the physical modes to fluctuate. The physical modes may be 
constructed by enforcing them to be invariant under the residual gauge transformations, 
(5a" = and 5hMr = (see appendix B). Thus, the physical fields we define are given in 
terms of the fluctuations, and classified as (from here on we set = k) 
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helicity two: E = gy^hyz, hyy - h^^ , 
helicity one: ^ = g^^ {ujhxy + khty)\ ay, 

and for helicity zero: 

1 i/c 9 k'^ n kuj n kw ^ 

$1 =4 - —af + -—a'l + — < + — 

(/) W(j) W(j) iO 

k'^fNwa^ k'^fw'^act) {aN' + 2Na') - Ir'^co^f {w^w' + k'^(t>') 
2r2w2 « + 4ruj'^w^ (/ + rf) 

2 i{-k'^ + w^) 3 
$2 =ai + ^ ^-at- 

UJW 

_ iA; 3 iw(t) irf{w^(l){aN' + 2Na') + 2N{P-w^)a(j)') 
" w""^ ~ + AuNwa{f + rf') ^^'^^ 

2k k'^fNa^ Ar^u^f -2ruj^f + k'^fa{aN' + 2Na') 

3*3 — ?x H ^tx 2 2 ^* o 2 / J- I Jh 

oj r'^LO'^ ircj^ (j + rj ) 

3 A; 3 w(j) I iojw 2 '^'^'t^ 

$4 =a^ H — a^ ^ —,a^ ^ —^a^ H 2 72^*^~ 

kfNw^a^cj) kf {fw^a^ {aN' + 2Na') + 2r^ {-u"^ + 0') 



with 



2r^uj (u;2 _ 02) 4^^ (^2 _ 02) ( J + rf,^ 



iy=gyyhyy, ^X = 9^"^ k^X , = ht , ^tX = Q'"'' hx ■ (3.4) 



First we look at the asymptotic behavior of the helicity zero physical fields in terms of 
the asymptotic form of the background (2.10) and the fluctuation fields (A. 5). The physical 
fields are chosen so that each one can be identified at the boundary with a fluctuation fleld. 
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dynamical fields 


constraints 


# physical modes 


helicity 2 


hyz ) hyy hzz 


none 


2 


helicity 1 


hty^ hxy] CLy 


hyr 


4 




htz ; hxz ) Oiz 


hzr 


4 


helicity 


httj hxxj hyy + hzz J hxtj CLf ' 


htn hxr: h^r^ CL^ 


4 



Table 1. Classifications of the fluctuations according to their transformation under the little group 
SO{2). The constraints are given by the equations of motion for the fields which are set to zero due 
the fixing of the gauge freedom: a° = and hrM = 0. The number of physical modes is obtained 
by the number of dynamical fields minus the number of constraints. Due to SO{2) invariance the 
fields in the first and second line of the helicity one fields can be identified. 



or a combination of them. In fact, in this limit they asymptote to 

$i(a;,r) {aD^, 



<I>3(a;,r)^(e.)[;-(C,)' 



(3.5) 

' 



Note that this computation was done in the A; = limit, since this is the relevant limit for 
this paper. The resulting correlators from the helicity zero modes will be written in terms 
of this physical fields. 

3.2. Equations of Motion, On-shell Action and Correlators 

In the following we will focus on the response exclusively due to time dependent pertur- 
bations, i.e. k^^ = (w, 0,0,0). In this case in addition to the 50(2) symmetry, P\\ parity is 
conserved which allows us to decouple some of the physical modes in the different helicity 
blocks. In this section we obtain the retarded Green's functions G of the gauge theory 
corresponding to the stress-energy tensor T^*^ and the currents Ja , defined as two point 
functions, as in 



G'"'^f"'{k) = -i^ 


1 dtd^x e" 




emT^^-{t,x),TP'^{QM) 




1 dtd?x e" 




0(t)([j,^(t,f),jao,o)]), 


G^Zik) = -i^ 


1 dtd^x e" 




0(t)([r^^(t,f),j,^(o,o)]), 


Grik) = -i^ 


1 dtd^x e" 




0(t)([j,^(t,f),r'^^(o,o)]). 



(3.6) 



Here T^^'^ and J a are respectively the full stress-energy tensor and current, which include 
the equilibrium parts of Sec. 2, {T^^) and {Ja), as well as the corresponding dissipative 
parts which arise due to the introduction of fluctuations in our model. 

We use the methods developed in the context of gauge-string duality to extract these 
Green's functions. First we determine the on-shell action at the boundary of the asymp- 
totically AdS space from which we can easily read of the Green's functions using the recipe 
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described in the seminal paper [44] and its generalisation to the case of operator mixing 
(c.f. [45]). 

We refer the reader to [19] for the treatment of the helicity one and two modes. Here 
we present the analysis of the helicity zero fluctuations. 

3.2.1. Helicity zero modes 

The equations of motion corresponding to these fluctuations are very lengthy, therefore, to 
guarantee readability, we omit them here. They can be derived by expanding the action 
(2.1) up to second order in the fluctuations and varying it with respect to the corresponding 
fields. 

Due to the parity symmetry i^j in the A; = case the modes split into two blocks, one 
transforming oddly (block 1) the other evenly (block 2) under Py. 

Block 1 - Parity odd 

The first block is composed by the modes {al,a^,a^,^tx}- The contribution of these 
modes to the on-shell action is^ 



con shell _ 1 f d'^^ / t c ' ''^"^ Ml' 2 2'^ rO^f^Na 3 3/ 

'^'hei.o. bi.i j ^ \-^^t.^t. + — [a^a,+a^a, 



^hei.o, bi.i j -^—^ y^j^f.txc.tx -T -rut^t j - a^a^ 

3r^ ^ 9 r^a^ 



2/V 



itx'-"^itx {w'a] + <t>'al)^ 



(3.7) 



=r-bdy 



which is divergent as we send rbdy — ^ oo. The divergence can be cured via holographic 
renormalization, i.e. the addition of covariant boundary counterterms that cancel the 
divergences without affecting the equations of motion (see appendix A). To obtain the 
boundary action we plug the field expansions at rbdy into equation (3.7). Since we have 
four fields satisfying four second order differential equations and three constraints (coming 
from setting hxr, aj., to zero) we are left with a total of five (8 — 3 = 5) independent 
parameters at the boundary, {(«! )g , (o? )o ■, (ox)o ■> > (?ta;)o} (see (A. 6)). There is some 
freedom in choosing the undetermined coefficients, however the present choice is convenient 
for the later use of the gauge/gravity dictionary. We express the renormalized on-shell 
action at the boundary in terms of these coefficients^. 



2 



'!'hel.o,bl.i-^2y (27r)n^'-M2^ ^'-/^'^ 4 ^""^^0 

- 2"^o(6x)o + -jT^ K)o («?)o + Zp—f^ (°*)o ("')o (^-^^ 
-^rz^ («?)o («')o + ("')o K)5 + 2«''^? («')o (^<-)o 

^Here and in other similar expressions ahead, the products are to be understood as evaluated on opposite 
values of the frequency, as is natural for a Lagrangian written in Fourier space. For instance, ^t^, aj 
would actually be 5ta:(— Lj,r) aj(aj,r). 

■^AU fields in the following boundary action are dimensionless, i.e. we pulled out rn- Wherever the context 
allows, we are sloppy with the notation and do not give the dimensionless fields new names. 
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As we discuss in appendix B, there is a residual gauge freedom left, which has to be taken 
into account to obtain physically sensible observables. Using the gauge transformations 
given in (B.ll) for A; = and setting Kt = = Aq = 0, since they do not affect the fields 
discussed in this block, we obtain the unique linear combination (up to an overall scaling 
discussed in the previous paragraph) 



$4 



3 lojaj + (t)a\ - wcpCtx 
at + w -7, ^ 



(3.9) 



Following [33], we rewrite the boundary action (3.8) in terms of gauge-equivalent fields, 
which guarantees that our solutions are gauge invariant. The set of allowed transformations 
is parametrized by three coefficients. The gauge equivalents to the fields which solve the 
equations of motion and constraints are 



aQa^ — iojKq — 4>Kq — \ujwKx 



iwAg + 0Aj , 



2 

) 



(3.10) 



ao6 



tx 



Note that we also took an overall multiplicative scaling factor into account. This can be in- 
cluded because different solutions of the equations of motion are related by a rescaling of the 
fields. These expressions give a relation, parametrized by four coefficients {oq, Aq, Aq, K^}, 
between different sets of classical solutions which are equivalent. 

In order to compute the two-point functions as derivatives of the classical action, we 
will follow the directions given in [33], which instructs us to prescribe the value of the 
perturbations at the boundary, respectively defined as /3tx}, in terms of the 

gauge-equivalent quantities defined in (3.10). Those are 



/3t = "0 iaj)[ 



Ptx 



- ao (a^jQ 



ao (^i 



tx)Q 



iwAn + 



(3.11) 



The four coefficients of the gauge transformation can be chosen so that the fields asymptote 
to these vales. Thus, we are effectively fixing the gauge, because the gauge freedom is 
"absorbed" in the freedom of choosing the boundary values. Then, we rewrite the boundary 
action (3.8) in terms of the /3i and obtain 



on-shell 
hel.O, bl.l 



d^A; 



(Pl\ 

Pi 

\Ptxj 



(3.12) 



where ^(i) is the Green's function matrix of this block, which relates the response of the 
system to field fluctuations a\, af, and htx- Note that, following our convention, the 
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fields on the left row vector are evaluated in —oj and the fields on the right column vector 
are in u. Next, by taking derivatives d'^ /df3*{—uj)df3{uj) of the action above we obtain 



V(r*-)(a;)y 



/r c-on-shcll \ 
"'-'hclicity ^ 

X con-shell 
"'^hclicity 

r con-shcll 
"'-'hclicity 

'5(*4)S(-^) 

r con-shell 
"'-^helicity 



/Gilo 



t,x I 



^1, 



G^lico) G^^ico) G'^{u) G^iu) 
Gtiiu:) Gliiu:) Gf^u) 
G^^liuj) G*^3(w) G*^'*^(a;)y 



('f4)^(^) 



which explicitly written in terms of the field theory expectation values is 



^2 


-L^2 




iuj 




















\ 







-iJi) 








■{Tu)/ 



(3.13) 



where we already included a factor of 2 coming from the prescription developed in [44] 
for real-time correlators. Note that using the prescription above we automatically get the 
correlator which includes the physical field <&4 instead of as it is pointed out in [33] . 

The matrix is completely determined by the background solution near the boundary, 
except for one entry, the one corresponding to the two-point correlator of ^4, which in 
terms of the parity odd helicity zero modes reads 



~lX,X I 



1 2 2^ 2a2 
-a uj H r 

2 (r.'i\^ 



(^i)o 



4) J + 



u a 



+ iu {al)\ 



(3.14) 



Rewriting this correlator in terms of the physical field we obtain 



^3,3 



-a 



Mr) 



<I>4(r) 



+ counter terms 



(3.15) 



='"bdy 



showing that all entries of ^(1) are gauge invariant. 

To compute this correlator we have to numerically integrate the equations of motion 
and constraint equations of this block. Since we choose infalling conditions at the horizon, 
we fix four of the eight independent coefficients. And out of the remaining four coefficients, 
three are fixed by the constraint equations, leaving us with one free parameter. This 
parameter corresponds to the overall scaling of the physical field and is related to qq 
in (3.10) and (3.11). Since the correlators are defined by ratios of the boundary values, 
this parameter is scaled out and we can just set it to one. From the solution of the numerical 
integration, we read off the boundary values of the fields and plug them into (3.14) to obtain 
the Green's functions. The results are presented in section 4, together with a qualitative 
analysis of the thermoelectric effect associated to these correlators. 
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Block 2 - Parity even 

The second block is composed by the modes {a^^a^^a^^^^t^iy^^x} , which combine to 
form three physical fields. The combinations we chose were defined in (3.3), and in this 
section we are taking /c = 0, in which case they reduce to 

^1 fw' , 



2 \w 3 iwcf) irfw{2(j)Na' + (l)N'a -2(l)'Na) 
+ " + AuNa if + rf) 



2r/' - / 
f + rf 

The contribution of this second block of helicity zero modes to the on-shell action is 
™heii _1 /■ /^^^^. ra'^fNa / . , r^g^ 

'-'hol.O, bl.2 —^2 J (^27r)^ 1 4 2 "r "x"!: / 2cT * 



=»'bdy 

(3.17) 



which again is divergent. The renormalized on-shell action is derived and presented in the 
appendix, see (A. 14). 

Since we have six fields determined by second order differential equations and three 
constraints, we end up with twelve (12 — 3 = 9) undetermined coefficients of the boundary 
expansion, in terms of which the expression above is written. They are 

Notice that six of them (the ones with subscript 0) coincide with the boundary values of 
the fields. The other three are higher-order coefficients. They are undetermined since the 
boundary expansion does not know about the boundary conditions set on the horizon (i. e. 
that they must satisfy an infalling condition at the horizon). Actually when integrating 
the equations these coefficients are fixed by the boundary values we choose at the horizon. 
However, how the expansion coefficients at the boundary depend on the coefficients at the 
horizon cannot be addressed analytically, since in the bulk we can only solve the equations 
of motion numerically and this dependence is precisely determined by the behaviour of the 
fields in the bulk. 
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As in the previous case, the Green's functions cannot be extracted directly from (3.17) 
because there is a residual gauge freedom left under which the fluctuation fields are not 
invariant. To fix the gauge freedom, we can apply again the formalism used before to derive 
the gauge-equivalent solutions. In this case, we have to look for the restricted set of gauge 
transformations and rescalings that keeps unaffected the perturbations of the first block. 
This set is parametrized by six coefficients {a^, Aq, Kt, Kr}, with i = 1,2,3, and gives the 
gauge-equivalents of a solution, which are 

aoaf - iwA^ + iuj(f>Kt + (\/iV0' - u^{4>A - C^)) Kr , 



— )• aoo^ — wAq — iujwCfpKr 
al. — )• aoai + VNw' Kr , 



it aoit + 2i6^i^f + ( ^^ ^^'^ - 2lo'^A^ Kr , (3.19) 



2^N{1 + rf) 
rf 



Otoiy H Kr , 



2VN{f-2rf 
c,x ao?x H 7 K, 



r ) 



rf 

where A, are defined in (B.7). We do not explicitly write the 3 independent scale factors 
out, rather we use a general scaling uq. As will be explained later, the 3 independent 
scale parameters are related to the freedom of choosing the value for 3 of the fields at 
the horizon. However, due to the complicated mixing of the fields in the bulk it is not 
known how this translates into the scaling at the boundary. Following the steps of [33], 
we would now proceed by prescribing the values of the perturbations at the boundary 
{15^,13^, 13]., f3t, I3y, l^x} by evaluating the asymptotic behavior of the gauge transformations, 
that we find is given by 

(3f={a^)l-iu;Al + iu;4Kt, 

/3t = (6)0 + 2Kr + 2iojKt , 
(3y = {Ql + 2Kr, 

Px = {ixt + . 

Since we do not know how the scale parameters enter the above equations, we have to alter 
our approach in deriving the action in terms of the physical fields. 

We parametrize the fluctuations in such a way that each physical mode asymptotes (see 
equation (3.5)) to the boundary value of a fluctuation fleld. For this reason, from here on 
we will work with 

ip{u,r) = ix{uj,r) + £,y{uj,r) , Crn{uj,r) = C^{u},r) - iy{uj,r) . (3.21) 
In addition, we perform a rotation of the a^., into 

at{^,r) = al{uj,r) +ial{uj,r) , a~ {uj , r) = al{uj , r) - i al{uj , r) . (3.22) 
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Accordingly, we rotate the corresponding physical fields into ^>-i- = <I>i it i ^2, so that 
their respective boundary values coincide with those of a^. This parametrization is more 
convenient, since the and a~ fields transform under the fundamental representation 
of the unbroken [7(1)3. That is, they behave in a similar fashion as electrically charged 
vector mesons do under the ?7(l)em- To make contact with the unbroken phase, we keep 
the parametrization also in the broken phase. Notice that these fields are conjugate of one 
another: (a^(w))* = o^(— w). 

Next, we will invert the definitions (3.16) and solve for the selected fluctuation fields 
(fi = {a^,^m} in terms of the corresponding physical fields $1 = {$±,^3}. The idea is 
to replace these three fields and write the on-shell action in terms of the physical fields of 
this block along with the remaining fluctuations ipi = {at,^t,^p}- This can be seen as a 
change to a more convenient basis, which guarantees that the resulting correlators are free 
of gauge ambiguity. 

We perform the replacement in (3.17) and in the corresponding counterterms (App. A. 2). 
In terms of the expansion coefficients of the physical fields at the boundary we obtain the 
on-shell action 



Qon-shell 
'-'hel.O, bl.2 



(2-)^ 



a' 

^2 L 



<I>4 (<!>_)? + (^-)o (^^ 



1 



1 



+ 



2fl — LO 

12a; 



64 



+ 



5/1 



+ 



19m^ 
96 



3J0 



UJ 

H — (jj 

b r 



32 



12(6)^' 



K)o 
2 iCtt 

fe2 



rrir 



16 



(^3)0 H^p)l 



2 i^t)'o 



(3.23) 



9iQl +12(Cp)S(e 



=^bdy 



The fields in this action are defined by equation (3.5) and (3.21) and below equation (3.22). 
This new action (including the <pi part), when written in terms of the block 2 perturbation 
modes, coincides exactly with what we have in (A. 14). 

The part involving the physical fields only can schematically be written as 



Qon-shell 
•-"heLO, bl.2 



[$*(-w, r)A{k, r)udr^j{uj, r) + $*(-a;, r)B{k, r)u^j{uj, r) 



where the derivatives dr^i evaluated at the boundary absorb the higher-order coefficients 
of the expansions (see (3.18)), in the same way that the $i(r;,) absorb the boundary values 
of the replaced fields. Of the matrices A, B; we only need to know their asymptotic values 
at the cutoff rbdy, which are given by 



A{oj, rbdy) 



12s 

-5« '^bdy 







12s 

-i«^bdy 







12 ' bdy , 



(3.24) 
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and 



\ 







24aj 



24a; 



a W^ 



(3.25) 



24<^ 



5 



33(^,'^bdyj 



with BssiuJ, r) = ^ [-4t^V2 + log (^) - 160/| - 19m^] . 

At this point we refer the reader to [45] for a prescription to calculate the Green's 
functions in systems where the operators mix. In appendix C, we discuss this prescription 
in more detail and show how it can be generalized to our case. 

The matrices A, B are the ones used to perform the numerical calculations described 
there, and the rest of the terms in (3.23) give directly the Green's functions. In fact the 
matrix of Green's functions of this block may be written as 



) \ 



(J5)(c.) 

{T^^ + Tyy){uj 



) J 



s Qon-shcll 
""-^hclicity \ 

'5(»?)o(— ) 

r Qoii-shcll 
"'-^helicity 

r con-shell 
"^helicity 

X oon-shell 
"^hclicity 

5(<I"3)o(-^) 

r con-shell 
"'-^helicity 

r con-shell 
"'-^hclicity 



^(2)( 



(^-)oM 

(^3)oM 



(3.26) 



where the entries are denoted by 

/ G^l,{u:) G^^iu) G^^iu) G^iu;) Gf{io) G\\u^)\ 



G%{u) Gl'^^u) G^4{u) Gl^{u) Gl^{u) (a;) 

Gf;3(a;) G?;+(a;) G!;_^(t^) ^^'"(l^) G?P(a;) G?*(a;) 

G™^(a;) G™^(w) G™5(tj) G™'™(tj) G™'P(w) G"^'*(u;) 

GP\{uj) G^liio) GP5(w) Gf'"(w) Gf'P(a;) GP'\oj) 

\G^\{oo) G%{u) G*!(a;) G*''^(a;) G*'P(w) G*'*(a;) / 

and we find that they are given by Q{2){'~^) = 



(3.27) 



/ 







2b (^f) 








\ 






Glfiiu) 


Gi'^:{u) 


G^"^(a;) 


~^z;rwi / 










G_";_"(a;) 


G?™(lj) 


8lu Wl / 









G™^(a;) 


G™?(w) 


G™''^(a;) 




-|((7^*)-2(r..)) 







M-<^ / fxy 




((7^t) - 2(r...)) 






V 









-|((7^t)-2(r..)) 




liTtt) J 
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Note that the index m is related to <^3 field, since ($3)0 = (^m)o- Furthermore the ^ indices 
are related to the and fields. Here all entries are functions of the background, except 
for the sector containing the physical fields and their couplings, which has to be computed 
numerically. As before, to do so we choose infalling boundary conditions at the horizon 
which leaves us with six of the former twelve free parameters. From the six left, three 
are fixed by the constraint equations. The value of the three remaining coefficients just 
scale the solutions and since we compute ratios of the boundary values it does not matter 
which value we choose for them. Note that a more detailed explanation of the numerical 
procedure we apply is described in appendix C. This part of Q(2) describes the dynamics of 
this block, and it is related holographically to several interesting properties of the superfiuid 
phase, as we describe in the next section. 



4. Transport Properties 

In this section we extract the transport properties of the holographic p-wave superfiuid 
from the correlation functions presented in the previous section. We split our analysis into 
distinct transport phenomena. 

4.1. Thermoelectric EflTect parallel to the Condensate 

We start by presenting the thermoelectric effect parallel to the condensate, i.e. we look at 
charge transport and temperature gradients in the x direction. This is related to the first 
block of helicity zero states we presented in section 3.2.1. Furthermore our results are in 
agreement with [33] for the non-backreacted case. 

The thermoelectric effect describes the simultaneous transport of charge and heat (or 
energy). This means that an electric field not only leads to a current, but also to a heat 
flux and, conversely, a temperature gradient leads to an electric current in addition to a 
heat flux. In holographic systems, this effect was already observed in s-wave superfluids 
(see e.g. [31, 37, 46]) or in the p-wave superfluid component transverse to the condensate 
[19]. However, in the case at hand, we have a slight complication due to a further coupling 
of the a\ and fields to the and to the htx metric component (see section 3.2.1). 

A straightforward calculation (see [37]) shows that VxT is related to the gu component 
of the metric through the change in the period of the Euclidean time. The change can be 
done in a way that 5gtt becomes pure gauge provided a complementary change is done for 
and gtx- It is customary to fix the gauge requesting 6gtt to vanish, and allowing for 
5A^ and 8gtx only. On the other hand. Ex receives an additional contribution from the 
vector potential, iuj (*I'4)o. The combined effect is so that we define the electric field and 
temperature gradient 



Ex = io; 
VxT 



(4.1) 



tx)o 



This modes source the charge current in direction of the condensate and the heat fiux 
Qx _ rptx _ ^jx ^ respectively. The relation of these currents to the corresponding electrical 
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field and temperature gradient defines the conductivity matrix 





Ex ' 

-(Vx.T)/r- 



(4.2) 



Comparing this matrix to the lower right corner of the one in (3.13), we can identify 
the electric, thermal and thermoelectric conductivities, which are related to the retarded 
Green's functions by 



--^3,3 ' 



rj,-XX 



i 



QXtX 



x,x 
3,3 



UJ 



(4.3) 



G 



tx,tx 



x,x 
3,3 



UJ 



((7^t)-2/x(J3*))+/iV 



2 XX 



The conductivity in direction of the condensate a^^ has been calculated numerically. The 
results are shown in figures 2 and 3 for a = 0.316 < Oc- The results for other values of a 
do not show any significant qualitative difference, therefore we do not show them in this 
paper. 

The rest of the matrix (3.13) shows the response of the system due to the aj, fluctu- 
ations. This is a manifestation of the fact that the equations of motion of the gauge field 
fluctuations are coupled. Therefore, if a temperature gradient excites one of these modes, 
the other two will respond, and their response is dictated by the coefficients of the matrix. 
In [33], the a^, of fiuctuation fields are interpreted as generating a rotation of the charge 
density in direction (J^) and (J2), however without changing its magnitude. 

The complete transport matrix of this block then reads 



(4) 
(J- 



( 



\- 



t.t 

^1,1 


t,t 
^1,2 


t,x 
0"l,3 


t.t 


t,t 
^2,2 


t,x 
0"2,3 


X,t 

^3,1 


x,t 
^^3,2 


^xx 


x,t 
/^^3,1 


x,t 

-^^3,2 





t,x\ 
-^^1,3 \ 

t,x 
-^^2.3 



rpj,XX 



I iuaj \ 



lujar 



V- 



Ex 

T 



(4.4) 



where each of the transport coefficients is simply related to the corresponding Green's 
function by o" = —iG/uj. We will now focus on the electric conductivity a^^, the others 
can be obtained from it. 

The fact that the longitudinal conductivity a^^ has a different behavior than that of the 
component transverse to the condensate a-^-^ (c.f. [19]) in the broken phase is an effect of 
the breaking of rotational symmetry. 

Let us discuss the similarities and differences between a^^ and o""'""'" (for a discussion of 
CT"'""'" see [19]). The curve of the real part of a^^ (fig. 2) shows the correct [47] asymptotic 
behaviour for large frequencies, i. e. the real part is proportional to the frequency for all 
temperatures. More precisely, for uj ^ T we have 



K52Re((T^^) 

2a2r 



UJ 



IT 



2'kT 



(4.5) 
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Figure 2. Real part of the conductivity Re(o'^^) over the frequency uj/{2-kT) for a = 0.316. The 
color coding is as foUows: blue T = 1.63Tc, red T = 0.98Tc, brown T = O.SSTc, green T = 0.78Tc, 
orange T = O.SOTc- There is a delta peak at strictly a; = 0, not noticeable in this figure, as dictated 
by the sum rule (the area below the curves has to be the same for any T). 

We expect this behavior on dimensional grounds, and as a consequence of the conformal 
symmetry in our system^. On the other hand, for decreasing frequencies we see that the 
conductivity decreases until nearly vanishing. This sharp decrease is a known feature of 
superconductors. It is present for all temperatures, not only for T < Tc- However, for 
smaller temperatures the decrease takes place at larger values of oo. It is by far not as 
sharp as in o""'""'", and furthermore there are some qualitative differences: The bump before 
decreasing is absent for a^^ - the asymptotic value for large frequencies is approached by 
the curves with smaller temperature from below, rather than from above, as opposed to 
the perpendicular case. Besides, up to numerical inaccuracy the conductivities do not seem 
to vanish for any frequency, for temperatures above 0.5Tc. In comparison, the transverse 
conductivity has a far stronger temperature suppression in the gapped region. However, 
below 0.5Tc the situation seems to change dramatically as is explained in the next section. 

The real part of o"^^, as opposed to the perpendicular case, increases again for small but 
finite frequencies and reaches a finite value in the w — >■ limit, as seen in the zoomed region 
of figure 2. This increase in the real part in the zero frequency limit is due to a quasinormal 
mode which moves up the imaginary axis in the complex frequency plane (see the blue arrow 
in figure 6) and seems to reach the origin w = at temperatures slightly above O.STc. The 
increase we see towards the a; — )• limit comes from the projection of the quasinormal mode 
onto the real frequency axis. Note that this bump increases with decreasing temperature. 
Unfortunately it is challenging to compute the exact temperature when the mode arrives at 
the origin, since we have to rely on numerical calculations. Nevertheless, for temperatures 

•^There is no lattice spacing which would spoil the high hequency behavior. 
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Figure 3. Imaginary part of the conductivity times the frequency, a;Im(cr^^) over the frequency 
uj/{2ttT) for a = 0.316. The color coding is as fohows: blue T = IMTc, red T = O.QSTc, brown 
T = 0.88Tc, green T = 0.78Tc, orange T = O.SOTc. The curves tend to a constant value as w ^ 0, 
which indicates the presence of a pole at the origin. This is related to the delta peak in the real 
part of a^^ . 

below 0.5Tc it appears that a pole is formed and at the same time the real part of the 
conductivity is more strongly suppressed at finite small frequencies in comparison to cases 
of temperatures above 0.5Tc (see the green and orange curve in zoomed region of figure 2). 
It seems that somewhere around 0.5Tc, due to the quasinormal mode at the origin, the 
conductivity behavior in the direction of the condensate changes. It would be interesting 
to understand this effect from a field theoretic point of view, we leave this for future work. 

Due to the pole in the imaginary part of the conductivity (see fig. 3) and the Kramers- 
Kronig relation [32] we know that at uj = the real part must have a delta peak. There are 
two main contributions to the prefactor of this delta peak, which change with temperature, 
and they come from the pole at the origin of the imaginary part, expressed as 



This is reminiscent of the perpendicular case (c.f. [19]). In fact, similarly the first contri- 
bution is a consequence of translational invariance at all temperatures, specifically for 
temperatures above Tc. The other contribution, A^, appears when temperatures decrease 
below Tc. This prefactor is expected to be connected to the superfiuid density, however it 
differs from the corresponding factor in the transverse case. 

The properties of the two components of the conductivities we state here are very similar 
to the ones found in the non-backreacted case (see [33]). Therefore corrections due to the 
backreaction seem to be rather small. 




(4.6) 
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Figure 4. These plots show the real and imaginary part of the correlators Gj'^ versus the reduced 
frequency uj/{2t:T) for a = 0.316 at different temperatures: blue T = O.gSTe, red T = 0.88Tc, 
brown T = 0.78Tc, purple T = 0.62Tc, green T = O.SOTc, orange T = 0.46Tc. 



4.2. Viscosities and Flavour Transport Coefficients 

The second block of coupled modes transforming as scalars under the SO (2) symmetry 
includes the fields a^, a^, a^, = g^*^htt, Cx = g^^hxx and = gy^hyy. Similarly to the 
first block, these fields form 3 physical modes, $i, $2 and $3 (see (3.3)). It turns out that 
it is more sensible to consider this fields in terms of = $1 it z$2, since they transform 
fundamentally under the ?7(1)3 in the unbroken phase. 

4.2.1. Piezoelectric effect 

The transport properties presented in this section show similarities to an effect known 
as piezoelectric effect^ found in crystals [39]. This effect describes the generation of an 
electric current due to the squeezing and/or elongation of a crystal, or the generation of 
a mechanical strain due to an external electric field. A coupling between a normal stress 
difference and (flavour) currents that resembles this effect is found in this block. 

In simple terms, there is an interaction analogous to the one of the first block, between 

'^A similar effect, the flexoelectric effect, is related to the helicity one modes, see [19] 
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(a) (b) 

Figure 5. These plots show the real and imaginary part of the correlator G^'^, or equivalently 
versus the reduced frequency ijj/{2ttT) for a = 0.316 at different temperatures: blue T = 0.98Tc, 
red T = 0.88Tc, brown T = 0.78Tc, purple T = 0.62Tc, green T = O.SOTc, orange T = OAQT^. 



The broken phase is characterized by a condensate (Jf). The fluctuate around this 
background value and the system reacts by working against this perturbations by changing 
the diagonal stress-energy tensor components (T^^), {T-^-^) and (T**). The converse case, 
where we fluctuate about equilibrium values of the stress-energy tensor and look at the 
response of the currents (Jf ) works in a similar fashion. Note that this is not the only 
response of the system to these fluctuations. However, in this section we are interested 
exactly in the coupling between different modes. In the field theory this may be related to 
an electric current being affected by, or generating, mechanical stress (Piezoelectric effect). 
Finally, the transport coefficients "measure" the strength of the response of the system, 
i.e. how do the expectation values change with respect to the original values when they are 
perturbed. 

In figures 4 and 5 we plot the real and imaginary part of G^Jj- and G^'^ over the reduced 
frequency oj/{2ttT), for several values of the temperature, or equivalently of the chemical 
potential /x. We find the symmetry relations 

G!;_^(w) = Gl'^^i-u)* , G+'!(a;) = G?;+(-t^)* , 

as we did in the study of helicity one modes [19]. This was expected from the fact that 
($+(a;))* = «'-(-^), = $+(-^) and (^-sM)* = ^3(-^)*. 

In figure 7 we plot the G™^, whose imaginary parts are identical to those of G^™, and 
whose real parts are similar, except for small frequencies compared to the temperature. 

Notice that many of the curves in figs. 4-7 show a pole at w = 0. To understand 
why, remember that the formation of (Jf ) selects a preferred direction in flavor space, 
spontaneously breaking the S0{3) and U{1)3 symmetries. As a consequence of this, the 

field becomes one of the three massless Goldstone modes arising from the spontaneous 
symmetry breaking. This common pole reflects the formation of this Goldstone mode, 
since it is included in the flelds ^± , which are involved in all of the correlators presented 
here. When plotting the quasinormal modes in the complex frequency plane we also see 
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Figure 6. This figure, taken from [34], shows the different quasinormal modes in the D3/D7 system 
in the complex frequency plane. Here tr = uj/{2ttT). The red and green curves show the modes of 
the fluctuations which correspond to $± in our setup and the blue curve corresponds to our $4. 
It is interesting to see that the backreaction and the bottom up approach we pursue in the system 
at hand behave in a very similar fashion as the D3/D7 probe setup. Note however, that there is 
one difference: due to the backreaction and consequently the rotational symmetry breaking in the 
superfluid phase, the decouple from the <I>±, contrary to what happens in the D3/D7 mode. 
Moreover, we only see the red and green modes in the $± sector and not in the sector (see [19] 
for a treatment of this modes). 



this pole at the origin, so that for T > Tc they asymptote to the origin of the frequency 
place a; = 0. In fact, although it is not apparent in figs. 5 and 7, the correlators vanish 
completely in the unbroken phase T > Tc, since in this limit the equations of motion of the 
bulk fields decouple. 

Another common feature of our results is the appearance of a rich structure for the 
correlators in the broken phase, including the formation of a bump located on the same 
value of the frequency for all of them. These bumps come from higher quasinormal exci- 
tations. With decreasing temperature, they move in the direction of smaller values of the 
negative imaginary parts and larger real parts of the frequencies. Therefore they become 
more accentuated with decreasing temperature. Nevertheless the quasinormal modes stay 
in the lower half complex frequency plane for all the temperatures we were able to check 
numerically. We will leave it for future work to investigate their behavior at finite spatial 
momentum and zero temperature. Note that their behavior is very similar to the one found 
in the D3/D7 model in [34] (see figure 6). Following [34], these bumps may be interpreted 
as bound states, e.g. mesons. However, since we do not have a precise knowledge of the 
the field theory side, this interpretation should be treated with care. We cannot say much 
more at this stage, without having the exact formulation of the hydrodynamics dual to this 
gravitational setup. 

The Green's functions G^' j_ (c.f. fig. 4) seem to have different asymptotic values, how- 
ever this is just a consequence of the small frequency range displayed here. Actually, 
they do asymptote to the same value for all temperatures in the limit of large frequencies. 
However this veils the interesting details at low frequency, therefore we do not show it 
here. Nevertheless, the large frequency limit is proportional to a;^, in agreement with the 
underlying CFT. 

Finally, note that the real parts of G^"^ and C"^ are not symmetric to each other. 



-23- 




27vT 



(a) 




27vT 



(c) 




27vT 



12000 
10000 
8000 
6000 
4000 
2000 







1 

1 






2 4 6 8 10 

(b) 



-6000 
-8000 









= J^Im(G-?) : 



(d) 



(e) 



10000 - 






5000 - 






-5000 - 






10000 - 









2 


4 6 8 10 

(f) 



Figure 7. These plots show the real and imaginary part of the correlators G^J™ versus the reduced 
frequency uj/{2t:T) for a = 0.316 at different temperatures: blue T = 0.98Tc, red T = 0.88Tc, 
brown T = 0.78Tc, purple T = 0.62Tc, green T = O.bOTc, orange T = OAGTc. We are not showing 
G"^^ because their imaginary parts are identical. Their real parts, however, show a different low 
frequency behaviour. 



In the latter one we see poles in the a; — )• limit. This poles are due to the fact that in 
the case we are dividing by the boundary value of which contains the (massless) 
Goldstone mode and therefore vanishes at w = 0, whereas in the case we divide 

by the boundary value of $3, whose quasinormal mode is not located at the origin. 
In the next subsection we look at G""'*", the Green's function generated by $3. 



4.2.2. Transport Coefficient associated to a Normal Stress Difference 

In the a; —7- limit the imaginary part of the two-point function of $3 asymptotes to a finite 
value different from zero, see fig. 8. We expect that this Green's function is related to the 
transport coefficient A defined in appendix D. In the presence of anisotropy, besides the two 
shear viscosities ijyz and rjxy, there are three other coefficients. But in a conformal fluid, 
two of them, (^^ and (y, vanish due to the tracelessness condition of the energy- momentum 
tensor in conformal theories. There remains an additional nonzero component, A, which 
may be related to the normal first stress difference as discussed below. 

If we consider perturbations around a background with unbroken SO (3) symmetry, i. e. 
with zero w{r), then the physical field ^>3 decouples from the other fields. Inserting the 
analytic solution of the AdS Reissner-Nordstrom black hole, 

ds^ = -N(r)dt'^ + ^T^dr^ + (dx^ + dy^ + dz^) , (4.8) 
N[r) ' 

its equation of motion can simply be written as 

-$3 + piV(r)$'3]' = 0, (4.9) 



N{r) 

where N{r) = r^ — ^^-l- ^g^f . This is the equation of motion of a minimally coupled scalar, 
thus we can apply the procedure developed in [24] to derive the value of the corresponding 
transport coefficient. The relevant part of the boundary action in the unbroken case for 
this mode is 



Qon-shell _ ' h ' ^ 



4 J {2^r 



^2$^(a;, r)a^$3(w, r) + B^z{k, r)^l{uj, r)^z{uj, r) 



(4.10) 



with i?33(a;, f) = -^ [— 4a;^r^ -|- 2a;^ log (^) — lOmg] and ttiq defined as in the AdS Reissner- 
Nordstrom solution. Using the result of [24] and the Kubo formula (D.8), we obtain for 
this particular case the viscosity coefficient 



s ~ 18 V 4^ 



(4.11) 



where s is the entropy density. In the dual field theory interpretation, this would be the 
transport coefficient associated to the first normal stress difference, which is the difference 
of the diagonal components of the stress tensor. This matter has been discussed in the 
condensed matter literature, in the context of granular jets [48, 49], polymeric solutions [50] 
and nematic crystals [51], but to our knowledge, it has not been studied using holography 
techniques. It is one of the progenitors of many interesting and unique flow-features, like 
the rod-climbing or Weissenberg effect. 

Liquid crystals and polymeric solutions are examples of states of matter which display 
a crystal-like orientation and, at the same time, are able to flow like a liquid, and they can 
be subject to deformations. The word nematic makes reference to a common phase of these 
materials, in which the molecules are self-aligned along a certain direction, even though 
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Figure 8. We plot Im(G""'™(w))/(w s) j^^^ over the temperature T/Tc for a = 0.032 (red), 
a = 0.224 (yellow), a = 0.316 (blue) and a = 0.447 (green). Note that a = 0.447 > and 
therefore the phase transition is first order leading to multiple values near the phase transition. All 
curves tend to l/(247r) at Tc, since in the unbroken phase $3 satisfies the equation of motion of a 
minimally coupled scalar. 

they do not acquire a fixed positional order. The formation of this phase spontaneously 
breaks isotropy. At the level of symmetries, these materials are very similar to what has 
been described in this paper. 

One of the distinctive phenomena observed in the flow of liquid crystal polymers in the 
nematic state is that of a negative first normal stress difference in shear flow. Generically, 
whenever a force is developed due to normal stresses in a compressible material, it will tend 
to push apart the two surfaces between which the material is compressed (e.g. normal radial 
compression on a cylinder is expected to lengthen its shape on the vertical direction). But 
in liquid crystalline solutions in shear flow, the opposite effect is observed. The Weissenberg 
effect is a manifestation of this property, which occurs when a spinning rod is placed in the 
material. Instead of being expelled outwards, the fluid is attracted towards the rod. 

This is translated into our setup as having a positive coefficient A. A positive strain 
difference h^x — h±± corresponds to a deformation that enlarges the longitudinal direction. 
One would expect the response to be the formation of a compressing normal stress in the 
perpendicular direction, which translates into a positive pressure difference {P± — Px), that 
is, a negative {Txx — T±±). We have derived a lower bound for this material close to zero in 
the unbroken phase, and the results indicate a positive coefficient x that is increased with 
io/T. Thus, the broken phase shows this quality of nematic crystals, and its effect becomes 
more pronounced for lower temperatures, or large frequencies. 
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5. Conclusion 

We have considered a holographic p-wave superfluid within SU{2) Einstein- Yang-Mills the- 
ory, in which the formation of the condensate spontaneously breaks the rotation symmetry, 
selecting a preferred direction x and keeping transverse isotropy in the y, z-plane. This 
remaining SO{2) symmetry group allows us to classify the perturbations about equilibrium 
according to their transformation properties, into three different helicity sectors. In this 
paper we have focused on the helicity zero states and studied their transport properties. 
Due to parity, the helicity zero sector splits into two blocks. In the first block, we 
obtain the thermoelectric conductivity in the direction parallel to the condensate, whereas 
the second block allows us to study the piezoelectric effect and transport properties related 
to the first normal stress difference. These are interesting new phenomena which are due 
to the anisotropy in our system. 

We see that the thermoelectric conductivity displays some differences with respect to 
the transverse case (i.e. the helicity one fluctuations), despite being qualitatively similar. 
In particular, the temperature suppression in the broken phase is much lighter in the 
presented longitudinal case for temperatures above around 0.5Tc. For temperatures below 
the temperature suppression of the real part a finite small frequencies increases dramatically 
and we see a pole at a; = 0. This is due to a quasinormal mode traveling up the imaginary 
axis of the complex frequency plane. At around 0.5Tc this mode arrives at the origin and 
stays there. 

On the other hand, in the parity even block we find a behavior reminiscent of the 
piezoelectric effect. Furthermore we see bumps in the correlators of this sector, which seem 
to be related to the generation of bound states. 

In the zero frequency limit, we find a non-zero value for the two-point function of the 
diagonal metric fiuctuations, which is related by a Kubo formula to a component of the 
viscosity tensor, denoted by A. Since A has the same dimensions as a shear viscosity we 
investigate its behavior by taking its ratio with respect to the entropy density s. We find 
that in the broken case A/s is temperature dependent, whereas in the unbroken phase, 
it acquires a small positive value for all temperatures T > Tc. The ratio X/s does not 
fall below the value we find in the unbroken case for all temperatures. The physical 
interpretation of this coefficient is the effect that an anisotropic strain gradient has over 
the normal stress difference, {T^x — T±±). We compare our result to what is found in the 
context of nematic crystal liquids [51]. 

We have determined the coefficients associated to these effects for generic values of the 
frequency and the temperature. Our results are valid as an effective description of the 
transport properties near the critical temperature Tc, where scale invariance is approached 
and simple models of AdS/CFT can be applied. 

For further progress, a detailed analysis of the hydrodynamics of anisotropic superfluids 
is desirable to give a further interpretation to our study. In addition, it would be interesting 
to perform an analysis at finite spatial momentum, which would allow us to investigate the 
dispersion relations of the normalizable modes [52] and to check if there are instabilities 
similar to the ones found in [53, 54]. 



-27- 



Acknowledgements 



We are grateful to Patrick Kerner and Amos Yarom for discussions. DF acknowledges fi- 
nancial support from 2009-SGR-168, MEG FPA2010-20807-C02-01, MEG FPA2010-20807- 
G02-02, GPAN GSD2007-00042 Gonsolider-Ingenio 2010, and ERG StG 306605 HoloLHG. 
DF is thankful to D. Mateos, J. Gasalderrey-Solana, A. Buchel, G. Manuel, J. M. Pons, J. 
Tarrfo and M.A. Valle for discussions, as well as to the Max Planck Institute for Physics in 
Munich and to Perimeter Institute for Theoretical Physics for hospitality during the initial 
and final stages of this work, respectively. This work was supported in part by The Cluster 
of Excellence for Fundamental Physics - Origin and Structure of the Universe. 

A. Holographic Renormalization 

The boundary part Set of the action (2.1) does not have any influence on the equations 
of motion, but it must ensure that the action is finite on-shell. It includes the Gibbons- 
Hawking boundary term and some additional terms that will constitute the counterterm 
action Set, needed to cancel out any divergences that may appear. Thus, the full action is 
written as 



24 



2 



i? - A - ^ F^^^F-Miv 



1 '.4 



+ ^ U'^X^K + Set, (A.l) 



where K is the trace of the extrinsic curvature. 

We will follow the lead of the references [55, 56] to perform the holographic renormal- 
ization and obtain the counterterm action. 

A.l. Asymptotic Behavior 

In this section we look at the behavior of the fluctuation fields {-F(r)} at the horizon and at 
the boundary. Eventually we will want to calculate real-time retarded Green's functions [44, 
57], therefore at the horizon, besides regularity^, we have to fulfill the incoming boundary 
condition. For this purpose the ansatz we plug in for the behavior of the fields near the 
horizon is 

nr)Ur,=elY.Pl'^H. (A.2) 

where Ch = r/r^ — 1, into the equations of motion of the fiuctuation fields. It turns out 
that, as expected, we obtain two possibilities for (3, namely 

with T being the temperature defined in equation (2.11). As said before, we choose the 
solution with the "— " sign which corresponds to the incoming boundary condition. The 
other solution represents the outgoing boundary condition. 



®The condition (j>{rH) = at tlie horizon guarantees regularity. Even with all fluctuations switched on, 
there is no need for any further constraint. 
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On the other hand, our ansatz at the boundary is similar to the one used for the 
background calculation in section 2. However, here we have to add a logarithmic term to 
get a consistent solution (c.f. [55]). Therefore we use 



Let us now use the above expansions for the helicity zero states (the expressions for the 
helicity one and two states can be found in [19]). In this case, the equations of motion 
for the fluctuation fields can be distributed into two blocks. In the first block, we have 5 
independent expansion coefficients at the boundary (8 free parameters from the 4 second 
order differential equations minus 3 free parameters due to the constraints). We choose 
them to be {^tx)Q , {o,t)o^ ('^i)o' (^^Oo ^^'^ {^x)i- the horizon, we already halved the 
independent parameters by choosing the incoming boundary condition. From the remaining 
4 parameters, we can get rid of 3 by using the constraint equations. Therefore, we are left 
with just one free parameter at the horizon. 

We can perform similar considerations for the second block. Here we have also 3 con- 
straints, but we are dealing with 6 fields, each with its corresponding second order differ- 
ential equation. Therefore at the boundary we have 12 — 3 = 9 independent parameters, 

namely {Cy)^ , (Cx)o , (6)o > (4)o' («^)o' ("t)o' ^^y'>2- horizon, 

as before, we already fixed 6 free parameters by choosing the incoming boundary condition. 
There are 6 — 3 = 3 free parameters that give a fully determined system. 

Now we will state the first few non- vanishing terms of the expansion at the boundary of 
the different fields, because we will need them later on to determine divergences in the on- 
shell action and to calculate the Green's functions. The explicit form of these expansions 
is 




(A.4) 



where = {rh/r)^ is the expansion parameter. 



{^yt + ^ 




[iivt-iiX^ log 6.) el + 0{el) 
\{i.t-{iyf^ log 6.) el + 0{el), 




1 



6 



48 




(A.5) 
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for the fields of the second block, and 



'Jo + 



4)0^ + 1^^ (at)o 



«t = la* 



)o + 



'0 



b) ' 



— 



^1 ia;(a3)J^u;J ^ 
— £6 i ^ 



^t.=(6.)o-«' {al)l4eb + 0{el); 



(A.6) 



for the first block. Note that = (j)Q, (f)\ and w\ are the expansion coefficients of cj){r) and 
w{r) at the boundary. 

We do not state the expansions at the horizon, since the explicit form is quite lengthy 
and does not provide additional information to equation (A. 2). 

A. 2. Counterterms 



By plugging the expansions (A.6) into (3.7) and (3.17), we obtain the non-renormalized 



on-shell action, S'on-shell = ;^ / , where the integrand Cr^ is written in terms of the 

free parameters of the previous expansions as 
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(A.7) 



This is evaluated at the boundary, where = (r/^/rbdy)^ = 0, so any higher order terms 
vanish. And we have changed to momentum space, so that in each product of expansion 
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parameters in this expression, the first one has always to be understood as evaluated on 
—uj and the second on a;, e.g. on the first term we have (a|)g(—a;) (0^)^(0;). Therefore, 
note that the order in which they are multiplied matters. 

The terms that have to be considered for the counterterms are the ones in (A. 7) with 
explicit Efo dependence, since those are the ones responsible for the divergences 



1 



^efe 8 
+ log efe 



2, ,2 



a ijj 



62 



a2(;x2-a;2) / ^2 



>b2 



2^6 




For the construction of the counterterms, first we need to define the induced metric 7^1/ on 
the r = rbdy plane, 



resulting in 



^x^^ dx'' 



-A^(rbdy)cT(rbdy)'dt2 + "^^^ 



f{rhdy) 



fbdy 



4dx2 + r2^y/(rbdy)'(dy2 + dz2). 



(A.8) 



(A.9) 



Note that the expansion of the determinant of the induced metric for r ^ 1 is divergent 
and is given by 



+ 



^b^ 



2M +8fe)S(6)S-2fe)S(6)S-M(6)S 



^b2 



48 r L 



288 



fe2 



+ '4(^t.)f + "4 (5(^^)0' + (6)0' + 4 {Ql fe)^ - 12 (e.)^ (6)^ - 6 fe)^ (6)1; 



(A.IO) 



We will use these divergences to cancel out the ones we find in the non-renormalized 
action, together with other counterterms that have to be considered. It is not necessary to 
rigorously derive the covariant counterterms here in this work. By looking at the ones that 
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B. Sahoo and H.-U. Yee calculated in [56], we get an idea of how they should look like; 
namely, some combinations of -R[7], Rfj,u[l\ and F^^^ (i.e. the Ricci scalar and Ricci tensor 
on the induced surface, and the field strength tensor on that surface). Possible covariant 
combinations of the three terms are -y/— 7, ^— 7i?[7], ^— 7i?[7]2, ^/—^R>^'^ [y]Rfiu [7] and 
yJ—^F'j^^F"'^^. The coefficients in front of them can be guessed by requiring the divergences 
to vanish in the complete action. Their expansions for r ^ 1 are 

,4 



4 



(^a;)o 



■ 2w^(a^)f - 2 + (^(ai)f + {al)^^^ + 8i/xw (ai)|; (a^)^ 

(A.ll) 

It can be checked that by adding the real space action 



5ct = y d^x (3 + \rYA + 



log eft 



(A.12) 

to the action S'on-sheii (2-1) we get a divergence-free theory (up to second order in the 
fluctuations) for rbdy ^ 1, i-e. also the real time Green's functions are divergence- free. 
The renormalized rbdy ^ 1 on-shell action of the helicity modes are 
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Since we have 6 fields determined by second order differential equations and 3 constraints, 
we end up with 12 — 3 = 9 undetermined coefficients of the boundary expansion, in terms 
of which the expression above is written. They are 

{(at)o ' (4)o ' ' («i)o ' ' (^*)o ' (^2^)o ' (?y)2 ' fe)o} • (A.15) 

B. Constructing the Gauge Invariant Fields 

The gauge group of the SU{2) Einstein- Yang-Mills theory can be a subject of formal 
studies, as outlined in [58]. It is shown that diffeomorphism- induced transformations of 
the metric functions and pure Yang-Mills transformations of the Yang-Mills fields ought 
not to be considered separately. On general grounds, we must look for the most general 
combination, which can be written as 

6Af = OtS* + ^'diAf - M'd,Y.\ 



M ' (B.l) 



SA^ = Ftij^ + ^ij^^ + ^"A^^ + ^^djA'i + diK" + e^^'M'AI ■ 

where the . . . indices denote the spatial coordinates {x,y,z,r} of our spacetime. The 
metric gij is the spatial metric and g^^ is its inverse, and the functions N and A/"* (called 
lapse and shift vector respectively) are defined as 

ds2 = gj^i^dx^^dx^ = -M'^dt^ + gij{dx' +M'dt){dx^ +AAMt). (B.2) 

A general infinitesimal gauge transformation acting on a perturbed solution is given in 
terms of the 8 descriptors {S*^, A"}. We define 

Qmn = 9MN + huN, . , 

AO. _ Aa . a ^ ' 

— + '^M- 

where gMN and ^4^^^ are the background fields of the hairy black hole solution that is 
considered in Sec. 2. This part of the fields is therefore fixed, and the fiuctuations Hmn 
and a'^j are our dynamical variables. Thus, for instance the variation of the fiuctuation 
field defined as = g^^hu is given by 5^t = g'^^Shu = g^^Sgu- 

Furthermore, since they are considered as perturbations, they will be of the same order 
as the parameters T,'^^ and A". This allows us to give simple expressions to their variations, 
which will be approximated to lowest order. 
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B.l. Residual Gauge Transformations 

In Sec. 3.1, we decided to choose a gauge where = and hMr = 0. This kind of gauge 
fixing is allowed as long as, for any given configuration, there exists a gauge transformation 
such that it makes these components vanish. Since there are 8 functions that categorize 
each possible transformation, in principle this is feasible. Here we will see that this is 
justified, however the gauge is not completely fixed by these choices. 

We begin by defining the background metric, which corresponds to the ansatz use in 
section 2, so it is of the form 

ds^ = QMNdx^^dy^ = -ciirfdt'^ + C2{rfdx^ + C3{rf{dy'^ + dz^) + airfdr^, (B.4) 

the only non-zero components of the background Yang-Mills field are = (/){r) and 
a], = w{r), and we will be working in momentum space, i.e. 

S^^(t,x,r)= / d^x e^'=^^"S^^(tJ,/fc,r), 
A«(t,x,r) = / d^x e^^^^"A''(tJ,fc,r), 

where = {uj,k, 0,0), since in the case we are studying the rotational symmetry 50(2) 
is preserved so that the fluctuations can be classified. 

With these assumptions, we look at the variations of the h^jr components of the metric 
and the a" components of the Yang-Mills field, under an infinitesimal gauge transforma- 
tion (B.l) acting on a perturbed background solution. To first order, these are 



Shtr 


= -iwc4^S'' + Cl'S* - CldrT.\ 


(B.6a) 


Shxr 


= ikc4^Tr + C2^drT.'', 


(B.6b) 


S hyf 




(B.6c) 


6hzr 


= C3'a,s^ 


(B.6d) 




= 2C4 (C4'S^ + C49^S^) , 


(B.6e) 


Sal 




(B.6f) 


5a^ 


= drh\ 


(B.6g) 


5al 


= dr(l) + drh?. 


(B.6h) 



It is easy to convince oneself that by choosing carefully the Y^m and A"' functions, one could 
make the hMr and a'^ vanish. Now the residual gauge freedom would correspond to any 
further transformation that, while keeping these components null, changes the rest of the 
dynamical variables. We will find the most general form of a residual gauge transformation. 



-34- 



The solutions to ShMr = 0, (5a^ = can be written in terms of 8 constants {Km, Ag} as 

C4 



T,^{uj, k, r 

A^(a;, k,r 
A'^{uj,k,r 
A^{uj,k,r 



—KfCi — iujKrCiA, with A 
- ikKrB, with B = [ dr 

Ky, 

Kr 
C4 

ikK^Cyj + Ag, with Cyj = 



dr- 



C4 



(B.7) 



dr 



C2^ 



-Kt(i)-iujKr{(t)A-C^) + kl, with 



dr 



The physics ought to be invariant under any gauge transformation. Therefore, those dy- 
namical fields affected by these residual gauge transformations must be unphysical. Those 
linear combinations with the property of being invariant constitute the physical fields. 

B.2. The Physical Fields 

The helicity two fluctuations, H = g^^hyz and hyy — h^z are already invariant, that is, 

0, 



6E = gyydh. 

5{hyy 

therefore they are already physical modes. The helicity one fluctuations transform as 



"yz 

hzz) = 0; 



5h. 



i/ccs Ky, 



xy 

5hty = -iujcs^Ky, 
6al = 0; 



(B.9) 



so that the are physical, and the invariant combination of the other two gives the 
physical mode ^' = gyy{ujhxy + khty). Note that the same applies to the z components, 
which behave exactly the same as the y components. 

Now, for the helicity zero flelds^ ^tx, ^t, Cx, ^y, a" '^u 'we arrange any possible phys- 
ical mode $ as a linear combination of them given by some r-dependent coefficients r^, so 
that its invariance translates into 



(B.IO) 



a=l 



^Where we had defined = g^^hyy, = g^'^'K^, = ff"/itt and ^tx = g'^'^htx in (3.4). 
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Each of the variations in this expression are given by 

6Ctx = -iwi^x + ik^Kt -iok(B + ^a) Kr, 
d^t = 2iu}Kt + (—- 2uj^a] Kr, 

\CIC4 ) 

dCx = 2ikK^ +(— + 2k^B] Kr, 
\c2C4 J 

OKy = J^r, 

C3C4 

da]. = ikKl + ikwKx + ^— + A;^ {wB - Cyj)^ Kr 

5a\ = — iwAj — — iLowKx — ojk {wB — Cw) K, 
60^ = ifcAg — wAq — iojwC^Kr, 
5al = -iwAo + Mo + ik(f)Cy,Kr, 
5al = i/cAg + wAI - ik(f)Kt + ujk - C^) Kr, 
,3 , ■ ,^T^ , f ^' , ,2 



(B.ll) 



6a^ = -iwAg + iu}(l)Kt + - w {(j)A - C^) j Kr- 

Plugging everything into equation (B.IO) results in 6 algebraic equations, due to the fact 
that the variation of the physical mode must vanish for any residual transformation, that 
is, for any Kt, K^, Kr, Aq, Aq, Aq. Thus, we can solve for 6 of the r„ coefficients in terms 
of the other four. The solution gives the most general gauge invariant combination and it 
turns out to be independent of the {A, B, C^, C^} functions. 

What we call the four physical fields, $j {i: 1, . . . , 4), are chosen as a set of independent 
fields that generate that invariant combination. There is more than one choice, but the 
one we have taken is 



ik k"^ ■:, kbj -J, kw 

tx~ 



wcp W(p OJ 



efNwa'^ k^fw^acf) (aN' + 2Na') - 2r'^uj'^f {w4>w' + k'^cj)') 



$2 =4 + 



2r'^uj'^ Aruj'^wcj) (/ + rf) 



LOW 



\k 3 iw(t)^ irf {wmaN' + 2Na') + 2N {k"^ -w^)a(l)') 

~ w""^ " + AuNwa{f + rf') ^^"^^^ 

2k k'^fNa'^ Ar^uj^f -2ruj^f + k'^fa{aN' + 2Na') 
UJ r^LO'^ 2roj'^ [j + rf) 

3^3 wcp 1 iuJW 2 '^'^4' 

$4 =ax + -at - TsOt - ,,,2 _^2«t + .,2_^2«^'- 



3 _^ 

,2^ 



kf^Nw^a'^cp kf {fw^a^iaN' + 2Na') + 2r^ {-oj^ + 



r6 + 



2r2w(w2 -02)''^ Aruiu^ -(p'^){f + rf') 
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C. Numerical Evaluation of Green's Functions 



Here we review and generalize the algorithm to evaluate Green's functions in cases when 
there is operator mixing [45]. The starting point of the algorithm would be a general 
bilinear bulk action for some fields $i(x^,r) given by 

S = j (fx dr Au{x,r) 9^$a + $i B^,{x,r) 8^,^, + Cu{x,r) (C.l) 

In principle, one could be considering a perturbed background solution, as in the prob- 
lem discussed in this paper, and it may be possible that there is some gauge freedom 
associated with those perturbation fields. In our case, that would be given by the trans- 
formation (B.ll). But of course, gauge symmetry implies that the only relevant fields are 
the gauge-invariant combinations of the perturbations (3.3). Therefore, the most sensible 
strategy would be to write the action in terms of these physical degrees of freedom <I>i, and 
proceed from there. 

C.l. Writing action in the correct basis 

Even though the action is constituted as a gauge-invariant itself, it may not be possible 
to express it in terms of the physical fields only. It depends on the number of fluctuation 
fields and the extent of the gauge freedom. Let's say, for instance, that after whatever 
gauge fixing, the perturbed background is described by N fields (pi{x) and we are left with 
a residual gauge freedom parametrized by M constants. Then, the set of the possible gauge 
invariant linear combinations of (^i(x) is generated by — M independent physical fields 
But the part of the action that is quadratic in perturbations may be of the form 



S= I (fx dr 



d^ipi aij{x, r) d^ifj + tpi r) d^ipj + (ft Cij{x, r) ifj , (C.2) 

with, assuming for simplicity dependence on r only, {2N + 1)A^ coefficients {aij^hij^Cij} 
(Note that a and c form symmetric matrices) . The only requirement upon this action is that 
it be invariant under any gauge transformation. This gives 2NM equations (one for every 
field or derivative of field, and for every transformation), from which some coefficients are 
determined, leaving (1 + 2N — 2M)N undetermined coefficients that one is free to choose. 
On the other hand, an action written using only physical modes is constructed using 
2(A^ - Mf + {N- M) = {l + 2N- 2M){N - M) coefficients. Thus, the freedom in writing 
a gauge-invariant action is always greater than what the <I>i allow for. 

We conclude that in general the action will not be expressible as in (C.l). Not with the 
$j being physical, gauge-invariant fields. The generalization of the algorithm consists in 
getting as close as possible to an expression of that kind, as we explain below. Our starting 
point, for now, will consist in taking the complete quadratic action (C.2) and forgetting 
about the gauge symmetry issues. Varying this action, one can obtain the equations of 
motion for the perturbation fields, integrate the Lagrangian by parts, insert the equations 
of motion to obtain the action evaluated on-shell and add the proper counterterms to 
cancel out any divergences. Finally, this expression can be transformed carefully into an 
integration in Fourier space (see [45] for a description of the procedure). 
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The physical fields obey a set of coupled equations of motion of their own, and each 
particular solution gives a vector of functions {<I>i(x,r)}, such that as we approach the 
boundary, it asymptotes to some boundary values. It is in general possible to normalize 
the physical modes and to parametrize the perturbations in such a way that each boundary 
value of a physical mode coincides with the boundary value of one of the fluctuations, 

$i(fc,r) >^,{k). (C.3) 

r— >oo 

So we can make an association one to one between the physical modes and n of the N 
fluctuation fields. At any other distance r, each physical mode will of course depend on 
the values all the other fluctuations that are involved in its definition. 

Then, the first of the instructions would be to normalize the physical modes and to 
choose the appropriate fluctuation modes in order to be able to make this association on 
the boundary. The second step is to invert the definitions of the physical modes and solve 
for the fluctuations ^pi{k) that enter in the association. The idea is to replace these fields^ by 
inserting that solution into the on-shell action. In doing so, one obtains a contribution that 
involves only the physical modes, another one with couplings between the physical modes 
and the remaining fluctuations, and finally some terms given in terms of these remaining 
fluctuations only. That is, Sq.s. = So.s.i + 5'o.s.2 where^ 

So.s.1 = j d'^k [«>i 2l(A;,r)„ dr^, + $i 5S(A;,r)„ , 

5*0.8.2 = j d'^k [^ia{k,r)ijdripj + ipib{k,r)ijdr^j + ^ic{k,r)ij(pj + (pidij{k,r)dripj + ^mj{k,r)(pj]^^, 

The associated ipi{x) flelds no longer enter in the action. Now let's assume that we cannot 
flnd an analytic solution to the n coupled equations of motion, which is expected except 
for some simple cases. Nevertheless, since this action is evaluated on the boundary r^, a 
possible analytic approach would be to solve for the equations of motion on the limit r ^ 
and obtain the asymptotic expansions of the fields. As shown in the expansions of Sec. A.l, 
the expanded solutions are usually not determined by the boundary values (^q only. There 
are also some undetermined coefficients ipp which can only be fixed by supplying initial 
conditions at a given point from which integration starts. Since these coefficients depend 
on the whole integration up to the boundary, they will not be solved for analytically. It is 
for this reason that some Green's functions can only be evaluated numerically. 

A convenient position to start the integration is the horizon of the bulk geometry, because 
the initial conditions can be made easily at that point by demanding incoming solutions. 
This condition is related to the fact that we will ultimately be calculating retarded Green's 
functions. For convenience, let us refrain here what has been stated elsewhere in the text: 
The condition at the horizon halves the number of degrees of freedom and from a basis of 
2N solutions (A^ fluctuations under second order differential equations), we end up with just 

^From this point on, indices I, J, . . . denote the n = N — M physical modes, while indices . . . denote 

the N — n fluctuation modes that have not been replaced. 
*We do not state it explicitly, but each term in the next actions includes the product of a field evaluated 

in k and another evaluated in —k. This is natural for a quadratic Lagrangian written in Fourier space. 
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N solutions. Furthermore, there are M constraints coming from the equations of motion 
of the gauge-fixed fields which reduce these solutions to the n = N — M degrees of freedom 
that manifest themselves through the physical solutions and can be found by requiring 
invariance under residual gauge transformations. However, the analytic expanded solution 
obtained around the boundary knows nothing about the incoming condition at the horizon. 
We have only the constraints, so accordingly we are dealing with a basis of 2N — M = N+n 
solutions - that is, N + n undetermined coefficients. N of them can be taken to be the 
boundary values fg. Therefore, the number of undetermined coefficients ipp is expected to 
be precisely n, the same as the number of physical modes. So a numerical integration of 
the equations of motion of the physical fields starting at the boundary is sufficient to fix 
them, since we are implicitly setting n initial conditions. 

The expanded on-shell action can be then divided into two terms, 5o.s. = 5*0.8.(1) + 'S'o.s.(ii) 
where 



So.s.{l) = J d"^^ 

5o.s.(ii) = / d'^x 



(C.4) 



provided that the expansions are arranged in such a way that the n ((/'Op do not cross 
with the N — n (99i)o, which in general can be done because the are fixed by integrating 
the equations of motion of the physical fields, which are obtained by varying the action 
with respect to the physical fields, and therefore the part of the action which contains the 
remaining fiuctuations (pi is irrelevant to them. 

In this expression, remember that ($1)0 = (vOq- Now, to obtain the Green's functions, 
the AdS/CFT prescription instructs us to take the functional derivative of the action with 
respect to the boundary values of the fields [44]. The Green's functions of the fiuctuations 
(fi can be easily extracted from (C.4) and read 

Gfj[k) = -A,„ Gfj{k) = -Ky (C.5) 

On the other hand, the Green's functions G^{k), associated to the boundary values of the 
physical modes cannot be extracted directly, nor can they be expressed in an analytic way. 
These are the most interesting Green's functions because they have physical meaning, as 
discussed in Sec. 4. We have computed them numerically following the method presented 
in the next section. 

C.2. Prescription for numerical solutions 

It is tempting to establish the identifications 50,8.(1) = 'S'o.s.i and 5'o.s.(ii) = >S'o.s.2- However, 
it is important to note that this is incorrect, because the equations of motion of the fluc- 
tuations are all coupled, and consequently the expansions of the ipi{k,r) may depend on 
boundary values of the replaced ifi{k,r). Thus, expanding 5o.s.2 it is possible to produce 
terms like the ones found in 50.8.(1). 

In order to deal with this, an effective action S^g can be constructed, using physical 
fields only, such that its expansion near the boundary reproduces exactly those undesirable 
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terms. This effective action is to be subtracted from 5*0. s. 2 to cancel them out, and at the 
same time added to ^o.s.i producing 

ScS + So.s.i = J A[$i(-A:,r)A(fc,r)„a,$j(fc,r) + $i(-/c,r)5(fc,r)„$j(fc,r)]^^^^ . (C.6) 

The matrices A, B are obtained analytically, but only their expression at the boundary is 
necessary. We refer to the action as "effective" because its contribution matches exactly 
5*0.8.(1)) but it is just an artifact - we do not derive any equation of motion from it. The 
fields that will be inserted are solutions of the equation of motions derived from the original 
action. 

These solutions are obtained by numerical integration, starting from some selected values 
of the horizon ($1)0. In fact, this set of values determines completely the coefficients of 
the expansion at the boundary. We may choose n linearly independent sets (<I>i)g^''^ = e[^\ 
in order to obtain n linearly independent sets of boundary values. In particular, a possible 
choice is 

eP) = (1,0,0,...), ep) = (0,1,0,...), e^^ = (..., 0, 0, 1) . (C.7) 

Alternate choices are possible. This is just the one we used because we got good numerical 
results (with less noise). A numerical integration can be performed for each set in order to 
obtain n independent solutions {<I>|''^(/c, r)}j extended in the bulk, which can be arranged 
in a matrix H{k, r), with entries 

H,,{k,r) = ^['\k,r). (C.8) 

Thus, the J*^ solution appears as the J^^ column. On the other hand, we know that when 
each physical field approaches the boundary, it asymptotes to the value of its associated 
perturbation, ((/7i)o {k). At any other distance or scale r, since the system of differential 
equations is coupled, they will in general evaluate to a linear combination of all the {(v5.t)o}j, 
so that the set of functions can be written as 

^,{Kr)=Fu{Kr){^,%{k). (C.9) 

In this way, all the dynamics of the fields is encoded in the solution matrix Fij[k, r), which 
has the nice property of becoming the identity at the boundary, Fij{k, ri,) = 5ij. 

Any complete set of n independent solutions to the equations of motion is enough to 
build the matrix F, because any solution (any one that satisfies the incoming condition 
at the horizon) can be written as a linear combination of them. In particular, the matrix 
F{k, r) must be linearly related to H{k, r) because each F^ column of F is composed by a 
set of solutions that asymptotes to ((/2j)o = for all J, except for {ipi)^ = 1. Since at the 
boundary, by definition, F is the identity, the linear relation must be given by 

F{k,r) = H{k,r) ■ H'^{k,rb) . (C.IO) 

This result enables us to calculate the solution matrix, which encodes the dynamics, from 
n numerically integrated solutions. Then, by inserting (C.9) into our on-shell action (C.6), 
we obtain 

W.0.1 = / d'^k M-k,r)J'uik,r)<^j{k,r)l^,^ , (C.ll) 
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where T = A drF + B F. But, since this is evaluated at the boundary, where 
the matrix F becomes the identity and the physical fields coincide with their associated 
fluctuations, we might as well write 

Ses+s.o.i = j d''k{ip,)l{-k)[A,^{k,n) drFUk,n) + B,,{k,n)]{^,)l{k). (C.12) 

The Green's functions can now be directly extracted from this expression using the AdS/CFT 
correspondence prescription, to give 

G^,{k) = -A{k, n)H'{k, n)H-\k, n) - B{k, n) . (C.13) 

Notice that this formula reproduces the well known result for the Green's function of a 
decoupled equation. Here, instead of the derivative of the field, there is the matrix of 
derivatives. And instead of taking the ratio with the boundary value of the field, a factor 
given by the inverse of the matrix of solutions is included. 

This completes (C.5), giving the way of calculating all the Green's functions of the 
problem. However, as opposed to the ones given in (C.5), these Green's functions are 
not determined by the background only, the solution to the equations of motion of the 
perturbed degrees of freedom enters through H. Their physical meaning is clearer and more 
important, since they correspond holographically to the correlators of the dual operators. 



D. General Remarks on Viscosity in Anisotropic Fluids 

The concept of viscosity is linked to the internal motion of a system that causes dissipation 
of energy [59]. In general, we may define a general dissipation function H, such that 
the dissipative forces that describe the internal motion are obtained from it as velocity 
derivatives. Typically, frictional forces are linear in velocities u^, which suggests that the 
general form of this function be quadratic in velocities. 

But, for an internal motion which describes a general translation or a general rotation, 
the dissipation is zero. Since it describes dissipative processes only, H ought to vanish for 
these configurations of velocities. Because of this argument, the function must depend on 
the velocities through the combination of gradients of velocities u^i, = ^ (V^Ui, + V^u^), 
rather than on the velocities themselves directly. Thus, the general form is given by the 
sum H = ^rj^'^'^PUf^jyUxp, where the coefficients rj'^'^^f define the viscosity tensor [38], whose 
symmetries are given by 

^fiuXp _ ^u^iXp _ ^fiupX _ ^Xppu ^ (D.l) 
The part of the stress tensor which is dissipative due to viscosity is defined by 

In the case of a fluid in the rest frame = 1, and in order to satisfy the Landau frame 
condition u^H^'^ = 0, the stress energy tensor (and the viscosity tensor, correspondingly) 
must have non-zero components only in the spatial directions i, j, . . . = {x, y, z}. In general, 
only 21 independent components of r/jj^/ appear in the expressions above. 
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For the particular case of an isotropic fluid, the tensor can be written using only 2 
independent components, which are usually parametrized by the shear viscosity rj and 
the bulk viscosity (, so that the dissipative part of the stress tensor can be expressed as 
ir^ = —2ri{u^^ — ^5''^u\) — (u\6^^ , which is a well-known result. 

For a transversely isotropic fluid, there are 5 independent components in the tensor r/*-''^^ 
Without loss of generality, we choose the symmetry axis to be the x-axis. The non-zero 
components are parametrized by 





= Cx-2A, 


j^yyyy 


^zzzz 


Cy- 




_ ^xxzz _ ^ ^ 


^yyzz 


= - 2 - 


- Vyz 


j^yzyz 


= f]yz ) 


^xyxy 


^xzxz 


Vxy ■ 



X 

n+^yz 



(D.3) 



So that the non-zero off-diagonal components of the stress tensor are 

n ^ — 'if]xyUxy , FT — '2'UxyUxz 1 ^ 

In this consideration we are including only the contribution to the stress tensor due to the 
dissipation via viscosity, and we find the terms in the constitutive equation which depend 
on the velocity of the normal fluid u^. But in general, there would also be terms depending 
on the derivatives of the Nambu-Goldstone boson fields = df^ip on the superfluid velocity 
and on the velocity of the director, which may contribute to the dissipative part of the stress 
tensor (the director is the vector pointing in the preferred direction). 

However, these terms do not contribute to the off-diagonal components of the energy- 
momentum tensor because (1) a shear viscosity due to the superfluid velocity leads to a 
non-positive divergence of the entropy current [59, 60], and (2) no rank two tensor can 
be formed out of degrees of freedom of the director if the gradients of the director vanish 
[61]. In our case, the second argument is fulfilled since the condensate is homogeneous and 
the fluctuations depend on time only. Even though these degrees of freedom will generate 
additional transport coefficients, they do not change the shear viscosities, so we can write 
Kubo formulae which give the shear viscosities in terms of the stress energy correlation 
functions. 

Following the steps given in [62, 63], we can derive the Kubo formula that relates the 
coefficient A. Let us consider a conformal fluid, so that Cx = Cy = ^ {can easily be 
shown using the tracelessness condition of the stress-energy tensor) and concentrate on the 
difference in diagonal terms of the dissipative part, which is 

On the rest frame of the fluid, substituting = (1, 0, 0, 0), we find that the stress difference 
is given by the Christoffel symbols, which for a perturbed, spatially homogeneous metric 
of the kind considered in this paper give 

n- _ uyy = -sx (f*, - f*,) = - g^'g^x dt (Cx - e.) , (D.6) 
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where = g^^hxx and = g^^hyy and in a metric with AdS asymptotics g Qxx = 1- By 
comparison with the expectation from the hnear response theory, this equation is directly 
related to the zero spatial momentum, low-frequency limit of the retarded Green's function 
of T^x — T^y, which is the operator associated precisely to Cm = — Cy, 

G"''"'{uj)= lim / dtd\e-^'='^^"'e(t)([r^^(t,x)-r%(t,f),r^^(o,o)-r%(o,o)]) . 

(D.7) 

The Kubo formula that gives the A viscosity is 

A = - lim Im G™'™(a;) . (D.8) 
Lu^o 3 uj 

References 

[1] J. Erdmcngcr, M. Haack, M. Kaminski, and A. Yarom, Fluid dynamics of R-charged black 
holes, JHEP 01 (2009) 055, [arXiv: 0809 . 2488]. 

[2] N. Banerjee, J. Bhattacharya, S. Bhattacharyya, S. Dutta, R. Loganayagam, et al.. 
Hydrodynamics from charged black branes, JHEP 1101 (2011) 094, [arXiv: 0809. 2596]. 

[3] R. Loganayagam, Anomalies and the Helicity of the Thermal State, arXiv: 1211 . 3850. 

[4] R. Loganayagam and P. Surowka, Anomaly /Transport in an Ideal Weyl gas, JHEP 1204 
(2012) 097, [arXiv: 1201.2812]. 

[5] D. T. Son and N. Yamamoto, Kinetic theory with Berry curvature from quantum field 
theories, arXiv : 1210 . 8158. 

[6] J. L. Manes and M. Valle, Parity violating gravitational response and anomalous constitutive 
relations, arXiv : 1211 . 0876. 

[7] D. T. Son and P. Surowka, Hydrodynamics with Triangle Anomalies, Phys. Rev. Lett. 103 
(2009) 191601, [arXiv: 0906. 5044]. 

[8] K. Jensen, R. Loganayagam, and A. Yarom, Thermodynamics, gravitational anomalies and 
cones, arXiv: 1207.5824. 

[9] S. Bhattacharyya, S. Jain, S. Minwalla, and T. Sharma, Constraints on Superfluid 
Hydrodynamics from Equilibrium Partition Functions, arXiv : 1206 . 6106. 

[10] K. Landsteiner, E. Megias, and F. Pena-Benitez, Anomalous Transport from Kubo Formulae, 
arXiv: 1207.5808. 

[11] T. Kalaydzhyan, Chiral superfluidity of the quark-gluon plasma, arXiv: 1208.0012. 

[12] L Gahramanov, T. Kalaydzhyan, and I. Kirsch, Anisotropic hydrodynamics, holography and 
the chiral magnetic effect, Phys. Rev. D85 (2012) 126013, [arXiv: 1203.4259]. 

[13] A. Gynther, A. Rebhan, and D. Steineder, Thermodynamics and phase diagram of 
anisotropic Chern-Simons deformed gauge theories, JHEP 1210 (2012) 012, 
[arXiv: 1207. 6283]. 

[14] D. Mateos and D. Trancanelli, Thermodynamics and Instabilities of a Strongly Coupled 
Anisotropic Plasma, JHEP 07 (2011) 054, [arXiv: 1106. 1637]. 

[15] M. Chernicoff, D. Fernandez, D. Mateos, and D. Trancanelli, Quarkonium dissociation by 
anisotropy, arXiv : 1208 . 2672. 



-43- 



[16] M. Ammon, J. Erdmenger, V. Grass, P. Kerner, and A. O'Bannon, On Holographic p-wave 
Superfluids with Back-reaction, Phys. Lett. B686 (2010) 192-198, [arXiv:0912.3515]. 

[17] M. Natsuume and M. Ohta, The Shear viscosity of holographic superfluids, Prog. Theor. Phys. 
124 (2010) 931-951, [arXiv: 1008.4142]. 

[18] J. Erdmenger, P. Kerner, and H. Zeller, Non-universal shear viscosity from Einstein gravity, 
Phys. Lett. B699 (2011) 301-304, [arXiv : 1011 . 5912]. 

[19] J. Erdmenger, P. Kerner, and H. Zeller, Transport in Anisotropic Superfluids: A Holographic 
Description, JHEP 1201 (2012) 059, [arXiv : 1110 . 0007]. 

[20] P. Basu and J.-H. Oh, Analytic Approaches to An-Isotropic Holographic Superfluids, 
arXiv: 1109.4592. 

[21] P. Basu, J. He, A. Mukherjee, and H.-H. Shieh, Hard-gapped Holographic Superconductors, 
Phys. Lett. B689 (2010) 45-50, [arXiv:0911.4999]. 

[22] P. Kovtun, D. T. Son, and A. O. Starinets, Viscosity in strongly interacting quantum field 
theories from black hole physics, Phys. Rev. Lett. 94 (2005) 111601, [hep-th/0405231]. 

[23] A. Buchel and J. T. Liu, Universality of the shear viscosity in supergravity, Phys. Rev. Lett. 
93 (2004) 090602, [liep-th/0311175]. 

[24] N. Iqbal and H. Liu, Universality of the hydrodynamic limit in AdS/CFT and the membrane 
paradigm, Phys. Rev. D79 (2009) 025023, [arXiv: 0809. 3808]. 

[25] A. Rebhan and D. Steineder, Violation of the Holographic Viscosity Bound in a Strongly 
Coupled Anisotropic Plasma, Phys. Rev. Lett. 108 (2012) 021601, [arXiv: 1110.6825]. 

[26] A. Buchel, R. C. Myers, and A. Sinha, Beyond eta/s = l/4pt, JHEP 03 (2009) 084, 
[arXiv:0812.2521]. 

[27] J. Erdmenger, P. Kerner, and S. Muller, Towards a Holographic Realization of Homes' Law, 
JHEP 1210 (2012) 021, [arXiv: 1206 . 5305]. 

[28] Y.-Y. Bu, J. Erdmenger, J. P. Shock, and M. Strydom, Magnetic field induced lattice ground 
states from holography, arXiv: 1210.6669. 

[29] A. Karch and E. Katz, Adding flavor to AdS/CFT, JHEP OQ (2002) 043, [hep-th/0205236]. 

[30] S. S. Gubser, Breaking an Abelian gauge symmetry near a black hole horizon, Phys. Rev. 
D78 (2008) 065034, [arXiv: 0801 .2977]. 

[31] S. A. HartnoU, C. P. Herzog, and G. T. Horowitz, Building a Holographic Superconductor, 
Phys. Rev. Lett. 101 (2008) 031601, [arXiv: 0803. 3295]. 

[32] S. A. HartnoU, C. P. Herzog, and G. T. Horowitz, Holographic Superconductors, JHEP 12 
(2008) 015, [arXiv: 0810. 1563]. 

[33] S. S. Gubser and S. S. Pufu, The gravity dual of a p-wave superconductor, JHEP 11 (2008) 
033, [arXiv: 0805. 2960]. 

[34] M. Ammon, J. Erdmenger, M. Kaminski, and P. Kerner, Superconductivity from 

gauge/gravity duality with flavor, Phys. Lett. B680 (2009) 516-520, [arXiv: 0810. 2316]. 

[35] P. Basu, J. He, A. Mukherjee, and H.-H. Shieh, Superconductivity from D3/D7: Holographic 
Pion Superfluid, JHEP 11 (2009) 070, [arXiv: 0810. 3970]. 

[36] M. Ammon, J. Erdmenger, M. Kaminski, and P. Kerner, Flavor Superconductivity from 
Gauge/Gravity Duality, JHEP 10 (2009) 067, [arXiv: 0903. 1864]. 



-44- 



[37] S. A. HartnoU, Lectures on holographic methods for condensed matter physics, Class. Quant. 
Grav. 26 (2009) 224002, [arXiv: 0903. 3246]. 

[38] L. D. Landau and E. M. Lifshitz, Course of Theoretical Physics, Volume 7, Theory of 
Elasticity. Pergamon Press, 1959. 

[39] P. de Gennes, The Physics of Liquid Crystals. Oxford University Press, 1974. 

[40] S. Kobayashi, D. Mateos, S. Matsuura, R. C. Myers, and R. M. Thomson, Holographic phase 
transitions at finite baryon density, J HEP 02 (2007) 016, [liep-tli/06 11099]. 

[41] V. Balasubramanian and P. Kraus, A stress tensor for anti-de Sitter gravity, Commun. 
Math. Phys. 208 (1999) 413-428, [hep-tli/9902121]. 

[42] S. de Haro, S. N. Solodukhin, and K. Skenderis, Holographic reconstruction of spacetime and 
renormalization in the ads/eft correspondence, Commun. Math. Phys. 217 (2001) 595-622, 
[hep-th/0002230]. 

[43] S. S. Gubser, F. D. Rocha, and A. Yarom, Fermion correlators in non-abelian holographic 
superconductors, JHEP 1011 (2010) 085, [arXiv: 1002.4416]. 

[44] D. T. Son and A. O. Starinets, Minkowski-space correlators in AdS/CFT correspondence: 
Recipe and applications, JHEP 09 (2002) 042, [liep-th/0205051]. 

[45] M. Kaminski, K. Landsteiner, J. Mas, J. P. Shock, and J. Tarrio, Holographic Operator 
Mixing and Quasinormal Modes on the Brane, JHEP 02 (2010) 021, [arXiv : 0911 . 3610]. 

[46] C. P. Herzog, Lectures on Holographic Superfluidity and Superconductivity, J. Phys. A A42 
(2009) 343001, [arXiv: 0904. 1975]. 

[47] R. C. Myers, A. O. Starinets, and R. M. Thomson, Holographic spectral functions and 
diffusion constants for fundamental matter, JHEP 11 (2007) 091, [arXiv: 0706 . 0162]. 

[48] T. G. Sano and H. Hayakawa, Simulation of granular jets: Is granular flow really a perfect 
fluid?, Phys. Rev. E86 (Oct., 2012) 041308, [arXiv: 1206.6932]. 

[49] M. Alam and S. Luding, First normal stress difference and crystallization in a dense sheared 
granular fluid, Physics of Fluids 15 (Aug., 2003) 2298-2312, [cond-mat/]. 

[50] N. Lyotard, W. L. Shew, L. Bocquet, and J.-F. Pinton, Polymer and surface roughness 
effects on the drag crisis for falling spheres, European Physical Journal B 60 (Dec., 2007) 
469-476, [arXiv: 0708. 4291]. 

[51] C. I. Mendoza, A. CoreUa-Madueo, and J. A. Reyes, Influence of Electric Fields and 
Boundary Conditions on the Flow Properties of Nematic- Filled Cells and Capillaries. 
InTech, 2012. 

[52] E. I. Buchbinder and A. Buchcl, The Fate of the Sound and Diffusion in Holographic 
Magnetic Field, Phys.Rev. D79 (2009) 046006, [arXiv:0811.4325]. 

[53] S. Nakamura, H. Ooguri, and C.-S. Park, Cravity Dual of Spatially Modulated Phase, 
Phys.Rev. D81 (2010) 044018, [arXiv: 0911 .0679]. 

[54] H. Ooguri and C.-S. Park, Spatially Modulated Phase in Holographic Quark-Cluon Plasma, 
Phys. Rev. Lett. 106 (2011) 061601, [arXiv: 1011.4144]. 

[55] K. Skenderis, Lecture notes on holographic renormalization. Class. Quant. Crav. 19 (2002) 
5849-5876, [hep-th/0209067]. 



-45- 



[56] B. Sahoo and H.-U. Yee, Electrified plasma in AdS/CFT correspondence, JHEP 11 (2010) 
095, [arXiv: 1004.3541]. 

[57] D. T. Son and M. A. Stephanov, QCD at finite isospin density, Phys. Rev. Lett. 86 (2001) 
592-595, [hep-ph/0005225]. 

[58] J. Pons, D. Salisbury, and L. Shepley, Gauge transformations in Einstein-Yang-Mills 
theories, J.Math.Phys. 41 (2000) 5557-5571, [gr-qc/9912086]. 

[59] L. D. Landau and E. M. Lifshitz, Course of Theoretical Physics, Volume 6, Fluid 
Mechanics. Pergamon Press, 1959. 

[60] C. Pujol and D. Davesne, Relativistic dissipative hydrodynamics with spontaneous symmetry 
breaking, Phys. Rev. C67 (2003) 014901, [liep-ph/0204355]. 

[61] F. M. Leslie, Some Constitutive Equations For Anisotropic Fluids, The Quarterly Journal of 
Mechanics and Applied Mathematics 19 (1966), no. 3 357-370. 

[62] D. T. Son and A. O. Starinets, Viscosity, Black Holes, and Quantum Field Theory, Ann. 
Rev. Nucl. Part. Sci. 57 (2007) 95-118, [arXiv: 0704. 0240]. 

[63] G. D. Moore and K. A. Sohrabi, Thermodynamical second-order hydrodynamic coefficients, 
JHEP 1211 (2012) 148, [arXiv: 1210.3340]. 



-46- 



